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Abstract 



We present ACS, NICMOS, and Keck AO-assisted photometry of 20 Type la supemovae (SNe la) from the 
HST Cluster Supernova Survey. The SNe la were discovered over the redshift interval 0.623 < z < 1.415. 
Fourteen of these SNe la pass our strict selection cuts and are used in combination with the world's sample of 
SNe la to derive the best current constraints on dark energy. Ten of our new SNe la are beyond redshift z = 1, 
thereby nearly doubling the statistical weight of HST-discovered SNe la beyond this redshift. Our detailed 
analysis corrects for the recently identified correlation between SN la luminosity and host galaxy mass and 
(— I ■ corrects the NICMOS zeropoint at the count rates appropriate for very distant SNe la. Adding these supernovae 

Q^' improves the best combined constraint on dark energy density, pde(^), at redshifts 1.0 < z < 1.6 by 18% 

(including systematic errors). For a flat ACDM universe, we find ~ 0.729^0 014 (68% CL including 
systematic errors). For a flat wCDM model, we measure a constant dark energy equation-of-state parameter 



. w = -1.013t° o?3 (68% CL). Curvature is constrained to 0.7% in the owCDM model and to - 2% in a 

■ model in which dark energy is allowed to vary with parameters wq and Wa- Tightening further the constraints 

on the time evolution of dark energy will require several improvements, including high-quality multi-passband 
photometry of a sample of several dozen z > 1 SNe la. We describe how such a sample could be efficiently 
obtained by targeting cluster fields with WFC3 on HST. 
^ ' The updated supernova Union2.1 compilation of 580 SNe is available at http://supernova.lbl.gov/Union 

, Subject headings: cosmology: general, supernovae: general, cosmological parameters, distance scale 
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1. INTRODUCTION 

More than a dozen years have passed since combined 
observations of nearby and distant Type la Supernovae 
(SNe la) demonstrated that the exp ansion of the Universe 
is accelerating at the current epoch ([Perlrnutter et al.l lI998 



Gamavich et aL|[T998t [Schmidt et al.llI998l: iRiess et af 



1998 



Perlmutter et al.l [19991) . While the significance of the result 
has been boosted with the inclusion of larger, better cali- 



Wood-Vasev etall [2007 



Kowalski et al. 200: 
2009bt iKessler et alj [2009; ,Guv et al.. ,201 



brated SN la data sets (|Kn op et al. 2003; Astie r et al. 200^ 

[Hicken et alj 
the cause of 



the acceleration remains unknown. Einstein's cosmologi- 
cal constant, for which w, the dark energy equation-of-state 
parameter, is exactly —1 and independent of time, is just 
one of several possible explanations that is consistent with 
the constraints from SNe la and the constraints from other 
probes, such as the Cosmic Microwave Background (CMB, 
iDunklev et al.l2009l:[Koma tsu et al.'201 V) and Bar yon Acous- 
tic O scillations (BAO, [ETsenstein et al. 2005; P ercival et alj 

SNe la constrain cosmological parameters through the com- 
parison of their apparent luminosities over a range of red- 
shifts. At the highest redshifts, z > 1, the Hubble Space Tele- 
scope (HST) has played and continues to play a key role, in 
discovering and confirming z > 1 SNe la dRiess et al.ll2004l 
[2007; Kuznetsova et al. 2008; Dawson et al J2009h and in pro- 
viding high-precision optical and near-IR lightcurves. While 
constant w can be constrained using z 0.5 SNe la, SNe la 
at z > 1 provide the necessary redshift baseline to constrain 
time-varying w and some astroph ysical systematics like inter- 
galactic dust ( iMenardet al.ll2010l) . 

Discovering and following distant SNe la with the HST 
requires substantial amounts of telescope time because the 
field-of-view is quite small compared to that of ground- 
based telescopes. Therefore, all HST SN la discovery 
programs have coupled the search for SNe la and their 
photometric follow-up with other scientific studies. The 
GOODS survey dPTckinsgjLe t al. 2003; Riess elal] 12004 
[2007t iKuznetsova eTSl [200g) is an example that provided 
a window to pr obe the high- z univer s e for studying galaxy 
evolution (e.g. iBeckwith et alj [20061; iBouwens etal] 120061; 
iBundv et alj|2005l) in addition to z > 1 SNe. Nevertheless, 
progress in building a large sample of z > 1 SNe la discov- 
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ered with HS T is slow. In the most recent SN la compila- 
tion (Union2. lAmanuUah et alj[2010() which consists of 557 
SNe after the lightcurve quality cuts, only 16 HST-discovered 
z > 1 SNe la were available to help constrain a time-evolving 
w. (Well measured ground-based z > 1 SNe account for an 
additional four) 

Targeting regions that are rich in potential SN la hosts, such 
as galaxy clusters, offers a more effective strategy for using 
HST for SN la studies. Some of the earliest SN la searches 
used this strategy when the field-of-view of ground-based im- 
ages were only a few arc minutes across. The first spectro- 
scopically confirmed high-redshift S N la was discovered in 
a gala xy cluster (SN1988U: z=0.31. [Norgaard-Nielsen et all 
1 1 9891) as was the fir st high-redshift SN la observ ed by HST 
(SN1996cl: z=0.83. [Perimutter et al.|[T998l [T999I) . However, 
the advent of large-format CCDs, CCD mosaics and imagers 
with wide fields of view quick ly led away frorn this a pproach 
in the late 80's and early 90's (iPerlmutter et al.lll991l) . 

Given the large increase in the number of very distant clus- 
ters that have been discovered over the last 10 years, the an- 
gular extent of these clusters, and the field-of-view that is 
available with HST, targeting galaxy clusters beyond z ~ 1 
is again an effective strategy. In this paper, we discuss results 
from our HST Cluster SN survey obtained using this strat- 
egy. In addition to increasing the yield of SNe la dis coveries 
per HST orbit by a factor of two (iDawson et al.ll20()9[) . we in- 
crease the yield of SNe la in early-type galaxies by a factor of 
approximately four (Meyers et al. 20 11). 

SNe la hosted by early-type galaxies offer several poten- 
tial advantages over SNe la found in a broader range of host 
types. Stars in early-type galaxies are considerably older and 
span a smaller mass range than stars in late-type hosts. This 
may lead to a more uniform progenitor population. Evidence 
for this can be seen in the distribution in light curve widths. 
SNe la in early-type host galaxies follo w a narrower distri- 
bution than SNe la in late-type galaxies (iHamuy et al.l 119961 
120001: IRiess etal J [19991: ISullivan et all l2006h . Interestinglv. 
the rel ationship bet\yeen S N la color and host galaxy type is 
weak ( ISullivan et alj|2010l) . Using data from the HST Cluster 
SN Survey, we confirm both of these relationships for z > 1 
SNe la in Meyers et al. (201 1). 

With the availability of larger, better-calibrated samples, ev- 
idence for a correlation between host galaxy properties and 
SN la luminosities after correcti ons for lightc u rve wid th and 
SN la color is now emerging. iHicken et al.l (l2009cl) found 
that SNe la in early-type galaxies (morphologically classi- 
fied as E and SO galaxies), are 0.14 ± 0.07 mag brighter 
after lightcurve shape and color corrections than SNe la in 
galaxies of later types. A relationship of roughly the same 
significance bet ween host galaxy mass' and Hubb l e resid uals 
was reported by Kell v et all ( l2010l) . ISumvan et al.] ( 120101) and 
[Lampeitl et al. (2010). Uncorrected, this relationship leads 
to a significant systematic error in determining cosmologi- 
cal parameters, as the fraction of SN la in galaxies with high 
specific star-formatio n rates increases with increasing redshift 
(ISullivan et alJl2010l) . We expect that the host mass correc- 
tion is a proxy for more profound physics behind the SN la 
explosion mechanism; SNe la in early-type galaxies may lead 
to a better understanding of this correlation, given that more 
accurate mass , metallicity and age can be assigned to early- 
type galaxies (iBruzual & CharlotlBOOSt iTremonti et"aDl2004l 



' Host galaxy metallicity, specific-star-formation rate, or age are also 
drivers, as these quantities are somewhat degenerate in current data. 
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lMarastonll2005h . 

An additional source of astrophysical uncertainty concerns 
the color correction that is applied to SN la luminosities. 
There appears to be at least two mechanisms for the redder- 
fainter relation: extinction from dust in the ISM, which must 
play a role at some level, and an intrinsic relation between 
color and luminosity due to the explosion itself or the sur- 
rounding environment. There is no reason to believe that 
the redder-fainter relationship should behave in the same way 
for both mechanisms at all redshifts, but the two effects have 
proven to be hard to disentangle. 

Early-type galaxies contain significantly less dust than late- 
type galaxies, so separating SNe la according to early and 
late types offers a way to study the intrinsic component and 
to perhaps estimate the relative contribution and importance 
of dust in a broader sample. An early-type only sample may 
also yield a Hubble diagram with smaller statistical errors. 
Early work, based on a few dozen SNe la without color cor- 
rection, suggested t hat SNe la in early-ty pe galaxies are bet- 
ter standard candles (ISullivan et al.ll20d3l) . The evidence from 
more recent works, which use larger samples and better data, 
revealed that SNe la exhibit intrinsic diversity in color, but 
support the original fin dings dSullivan et al. 2003) with lower 
statistical significance (ISullivan et al.ll2010l) . 

Per unit stellar mass, SNe la are far less common in pas- 
sive, e^l^^LZELSSlS^iSLS'^JiLSi^liSHE'ESiiS^^ galax- 
ies jMannucci et al.1 120051: ISulUvan et al.ll2006h . Depending 
on the way hosts are classified, about one in five SNe la at 
low redshift will be hosted by an early-type galaxy. At higher 
redshifts, the fraction is expected to decrease, due to a combi- 
nation of an increase in the amount of star formation and ob- 
servational selection biases. Galaxy clusters, which are rich in 
early-type galaxies, even up to z ^ 1.4, are an effective w ay 
of finding SNe la in early-type hosts (iDawson et al.ll20d9l) . 

This paper is one of a series of ten papers that report 
supernova results from the HST Cluster Supernova Survey 
(PI: Perlmutter, GO- 10496), a survey to dis cover and fol- 
low S Ne la of very distant clusters. Paper I (iDawson et al.l 
12009^ describes the su rvey strategy and discoveries. Paper 
II (Barbarv et all 120101) repor ts on the SN la rate in clusters. 
Paper III (iMeyers et al.ll201lT) addres ses the properties of the 
galaxies that host SNe la. Paper IV (iRipoche et al.ll201 ih in- 
troduces a new technique to calibrate the "zeropoint" of the 
NICMOS camera at low counts rates, which is critical for 
placing NICMOS-observed SNe la on the Hubble diagram. 
The current work. Paper V, reports the SN la lightcurves and 
cosmology from the HST Cluster SN Survey program. Paper 
VI (Barbary et al in prep) wi l l repor t on the volumetric field 
SN la rate. [Melbourne et al.l ( 120071) . one of several unnum- 
bered papers in this series, present a Keck A daptive Optics 
observation of a z = 1.31 SN la in H-band. iBarbarv et al.l 
( |200^) report the discover y of the extraordinary luminous 
supernova, SN SCP06F6. iMorokuma et al.l (120 1 Ol) presents 
the s pectroscopic follow -up observations for SN la candi- 
dates. iHsiao et akl (1201 ll) develop techniques to remove prob- 
lematic artifacts remaining after the standard STScI pipeline. 
A separate series of papers, ten to date, re ports on clus- 
ter studies from the HST Cluster SN Survey: Brodwin et al. 
(12010); Eise nhardtet al. (200 8); Jee et al. (2009); Hil ton et al.. 
J2OO7, 2009h: [Huang et al.1 ( 2009); Sant osinn (l2009h; 
[Sfi-azzullo etaLl ( 120 1 Ol) ; iRosati et al.. (.20091) ; and iJee et al.1 
(1201 II) . 

This paper is organized as follows. In Sj2] we describe the 
HST Cluster SN Survey, the search strategy and discuss SN la 



typing. In [J3] we describe the procedures we used to process 
data and present the SN la photometry. In [2] we update the 
Union2 sample by adding the new SNe la from this paper, 
and we use the revised compilation to constrain cosmological 
parameters in fjs] 

2. SN DISCOVERIES AND DATA 

The HST Cluster Supernova Survey targeted 25 high- 
redshift galaxy clusters in the redshift range 0.9 < z < 
1.5 with the ACS camera on HST. Cluste rs were selected 
from the IRAC Shallow Cluster Survey (lEisenhardt et al.l 
20081). the Red-Sequence Cluster Surveys (RCS and RCS-2, 
Gladders & Yee'2005VGilbank etalj20l3), the XMM Cluster 
Survey (Sahlen et al. 2009), the Palomar Distant Cluster Sur- 
vey (JPostman et al. 1996), the XMM-Newton Distant Cluster 
Project (iBohringer et al.ii2005l). an d the ROSAT Deep Cluster 
Survey rRDCS. lRosati et al.llI999l) . At the time we conducted 
our survey, the sample represented a significant fraction of the 
known z > 0.9 clusters. Here, we summarize the SNe discov- 
ered in our survey. 

2.1. SN Sample 

As described in IDawson et al.l (120091) . the survey produced 
a t otal of 39 like l y SNe during the active phase of the search. 
In IBarbarv et al.l (120 1 Ol) . types are determined for 29 of these 
candidates. (The remaining 10 do not have enough light 
curve information to determine type, since they lie outside of 
our fiducial search time window or our signal-to-noise cuts.) 
Twenty SNe are classified as SNe la, with confidence levels 
of secure, probable or plausible. A secure SN la is one that 
either has a spectrum that directly confirms it to be a SN la 
or one that satisfies two conditions: (1) it occurred in a host 
whose spectroscopic, photometric and morphological prop- 
erties are consistent with those of an early-type galaxy with 
no detectable signs of recent star formation, and (2) it has a 
lightcurve shape consistent with that of a SN la and incon- 
sistent with all other known SN types. A probable SN la is 
one that does not have a secure spectrum but satisfies one of 
the two non-spectroscopic conditions that are required for a 
secure classification. A plausible SN la is one that has an in- 
dicative lightcurve but we do not have enough data to rule out 
ot her types. Details of t he classification scheme can be found 
in IBarbarv et al.l (I20I0I). a nd de tails of the galaxy typing can 
be found in lMevers etaTI (1201 II) . 

Sixteen SNe are classified as either secure or probable. We 
use these SNe in the cosmological analysis. We include the 
photometry and lightcurves of an additional four plausible 
SNe la to illustrate the quality of the data and the potential 
for a similar sample with complete classification (and because 
additional host galaxy data may later bring one of these into 
the larger sample). Secure, probable and plausible SN la are 
listed in Table[T] together with their position, redshift and typ- 
ing. Postage stamp images of the SNe and host galaxies are 
shown in Figure [Tf 

We labeled each of our 25 clusters with a letter from 'A' to 
'Z' (excluding 'O' to avoid confusion with zero) and assigned 
supernova names as 'SCP'H-[discovery year]H- [discovered 
cluster]H-[SN ID]. The c luster IDs, coordinate s , and redshifts 
are fo u nd in tables i n IDawson et aH (120091) ; IBarbarv et al.l 
( I20 1 Ol) ; [Meyers et aTl (l20 Il|). T he cluster membership is dis- 
cussed in iMevers et al.l (1201 II) in detail and summarized in 
Table[T]along with host type information. 

Several SNe la deserve special mention: 
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Table 1 

Supernova from HST Cluster Supernova Survey 



SN name 


Nickname 




-^cluster 


RA (J2000) DEC (J2000) 


E(B-V)'' 


Confidence 




SNe Hosted by Cluster Early-Type Galaxies 


SCPOSDO" 


Frida 


1.014 


1.017 


02:21:42.066 -03:21:53.12 


0.025 


secure 


SCP06H5 


Emma 


1.231 


1.241 


14:34:30.140 +34:26:57.30 


0.019 


secure 


SCP06K0 


Tomo 


1.415 


1.414 


14:38:08.366 +34:14:18.08 


0.015 


secure 


SCP06K18 


Alexander 


1.411 


1.414 


14:38:10.665 +34:12:47.19 


0.014 


probable 


SCP06R12 


Jennie 


1.212 


1.215 


02:23:00.083 -04:36:03.05 


0.026 


secure 


SCP06U4" 


Julia 


1.050 


1.037 


23:45:29.430 -36:32:45.75 


0.014 


secure 


SNe Hosted in the Cluster 


SCP06C1" 


Midge 


0.98 


0.974 


12:29:33.013 +01:51:36.67 


0.019 


secure 


SCP06F12 


Caleb 


1.110 


1.110 


14:32:28.749 +33:32:10.05 


0.010 


probable 


SNe Hosted by Early-Type Non-Cluster Members 


SCP05D6 


Maggie 


1.315 


1.017 


02:21:46.484 -03:22:56.18 


0.025 


secure 


SCP06G4'' 


Shaya 


1.350 


1.259 


14:29:18.744 +34:38:37.39 


0.015 


secure 


SCP06A4 


Aki 


1.192 


1.457 


22:16:01.078 -17:37:22.10 


0.026 


probable 


SCP06C0 


Noa 


1.092 


0.974 


12:29:25.655 +01:50:56.59 


0.020 


secure 


SNe Hosted by Late Type Galaxies 


SCP06G3 


Brian 


0.962 


1.259 


14:29:28.430 +34:37:23.15 


0.015 


plausible 


SCP06H3" 


Elizabeth 


0.850 


1.241 


14:34:28.879 +34:27:26.62 


0.019 


secure 


SCP06N33 


Naima 


1.188 


1.026 


02:20:57.699 -03:33:23.98 


0.023 


probable 


SCP05P1 


Gabe 


0.926 


1.1 


03:37:50.352 -28:43:02.67 


0.011 


plausible 


SCPOSPg" 


Lauren 


0.821 


1.1 


03:37:44.513 -28:43:54.58 


0.011 


secure 


SCP06X26 


Joe 


1.440 


1.101 


09:10:37.888 +54:22:29.06 


0.019 


plausible 


scpoezs" 


Adrian 


0.623 


1.390 


22:35:24.967 -25:57:09.61 


0.021 


secure 



SNe with No Definitive Redshift Measurement 



SCP06E12 Ashley ■■■ 1.026 14:15:08.141 +36:12:42.93 0.009 plausible 

^ Spectroscopically confirmed as a SNe la 

•> Redshift from SNe la or host galaxy IMorokuma et al.l2010i Barbary et al. 2010, Meyers et al. 201 1) 
' Redshift from cluster (Meyers et al. 2011, references therein) 
Galactic Extinction from lSchlegel et alj 419981) 

SN SCP06C1, a supernova discovered in 2006 in cluster C 
(XMMU J2205. 8-0159) and numbered as '1' among other 
transient candidates, could not be clearly associated with a 
galaxy in the cluster but was spectroscopically confirmed as 
a SN la at the cluster redshift. It might be an example of a 
SN la that comes from a pro genitor in the intracluster stellar 
population jSand etalj|2010l) SN S CP06C1 is discussed in 
greater detail in [Barbary et al.l (120101 ). As there was a bright 
background galaxy near the position of the supernova, had 
we not obtained a spectrum of the supernova, we would have 
misidentified the redshift. We note that experiments that do 
not obtain spectroscopy of their supernovae and instead as- 
sume the nearest visible galaxy is the host will have to factor 
cases like SN SCP06C1 into their analyses. 

SN SCP05D6, SN SCP06G4, and SN SCP06N33 occur 
behind the cluster, and are therefore gravitationally lensed 
by the cluster In most surveys, the lensing for SNe la 
hosted by field galaxies averages to nearly zero. In this 
survey, we target regions with larger than average magni- 
fication, so we must make the correction for SNe behind 
the clusters. To estimate the amount of magnification, we 
use the virial mass, M 200, from our weak lensing measure- 
ments (iJee et alJ lIoTlh . We assume a spherica l Navarro- 
Frenk- White profile (NFW; iNavarro et all Il997b with the 
concen tration parameters determined by the M200 — c rela- 
tion m iDuffv et all (I2OO8I) . The lensing magnifications for 
SN SCP05D6, SN SCP06G4 and SN SCP06N33 are esti- 



mated to be imitoml, 1-O15to:oo4 and 1.066^!];°^^ respec- 
tively. The magnification of SN SCP06N33 is larger than the 
others because the cluster is massive, and the host galaxy is 
located at approximately half the distance from cluster center 
as the others. We apply these corrections when using these 
SNe la in the cosmological fits and propagate the uncertain- 
ties accordingly. 

SN SCP06C0 is a more interesting case. The host is also be- 
hind the lensing cluster, which means that SN SCP06C0 will 
be lensed as well. Following the methodology used to cor- 
rect SN SCP06G4, SCP SN05D6 and SCP SN06N33, we find 
the magnification of SN SCP06C0 to be imotoml- Upon 
closer inspection of the host in the stacked ACS data and in 
more recent WFC3 data, a second much fainter object, pro- 
jected only 0'.'6 from the center of the host and about 0'.'2 
from SN SCP06C0, was detected. The object could be a satel- 
lite of the host, or it could be an unrelated galaxy along the 
line of sight. A spectrum of SN SCP06C0 was taken when it 
was about 6 days after maximum light; however, the signal- 
to-noise ratio was insufficient to allow a clear detection of 
light from the supernova given the observing conditions and 
the relative brightnesses of the hosts and the supernova. The 
classification and redshift therefore rely on correc tly assign- 
ing the supernova to the brighter galaxy. In ii l3.2.3l we model 
the surface brightness distribution of both galaxies. At the 
location of the supernova, the surface brightness of the large 
galaxy is four times greater than the surface brightness of the 
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Figure 1. Composite color (4775 and zgso) images of 20 SNe la from the HST Cluster Supernova Survey. Each SN la is shown in a box of 3.2" X 3.3" (North 
up and East left). Note the redshift of SCP06E12 is uncertain, and we use the cluster redshift as a guide. 
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small galaxy. Since [Oil] was not detected in spectra that were 
taken after the supernova had faded from view, neither galaxy 
is actively forming stars. The relative supernova rate is there- 
fore directly related to relative surface mass density of the the 
two galaxies. Using the surface brightness as a proxy for the 
surface mass density, we therefore assign SN SCP06C0 to the 
larger of the two galaxies with ^ 80% confidence and include 
this SNe la in the cosmological fits. 

3. PHOTOMETRY 

In this section, we describe the steps that were used to pro- 
cess the ACS and NICMOS data after they had been processed 
with the standard STScI pipelines. For the ACS data, we re- 
moved the spatially variable background from the pipeline 
processed data and applied charge transfer efficiency (CTE) 
and red-halo scattering corrections to the extracted fluxes. 
For the NICMOS data, we processed the data to compensate 
for amplifier offsets, bright Earth persistence, contamination 
from the passage of the telescope through the South Atlantic 
Anomaly (SAA), residual amplifier glow and fringing, and 
applied a wavelength-dependent non-linearity correction. A 
more detailed description of the individual steps now follows. 

3.1. ACS Processing and Photometry 

In general, the search consisted of four zs5o exposures 
per epoch and a single i775-band exposure. AU expo- 
sures are geometrically corrected (and multiple exposures 
are stacked) using M ultiDrizzle (Fruchter & Hook 2002 
iKoekemoer et al.l 120021) . The photometry is performed on 
the stacked zgso-band image and the single zyys-band driz- 
zled image. In addition to the data that w ere taken dur- 
ing th e search, five clust ers, CL1604-H43 04 (iPos tman et al.l 
l200l. RDCS0910-H 54 (iMei et all Hw^. RDC S0848-H44 
(IPostman et alJ2005h . RDCS 1252-29 (iBlakesleeet al.l2003l) 
and XMMU2235.3-2557 ( iJee et al.l |2009|) were observed 
with ACS prior to our program (PID9290 and PID9919). We 
have included all these data and processed them and the search 
data in a uniform way. In total, 1006 ACS exposures were 
processed. 

After processing the ACS data with the standard STScI 
pipeline using the most up-to-date calibration files, we re- 
moved the spatially variable background by masking all ob- 
jects and artifacts, and subtracting a heavily smoothed, me- 
dian version of what remained. Given the time baseline of our 
observations, guide stars changed between epochs. We must 
therefore use objects in the images - typically 20 to 30 objects 
per image - to tie the relative astrometry between epochs. The 
size of the residuals was typically 0.2 pixels, which is larger 
than one usually expects from point sources with good signal- 
to-noise ratios. We attribute this to temporal and spatial un- 
certainties in the distortion correction that are applied to im- 
ages that are both temporally separ ated by many mo nths and 
rotated with respect to one another (I Andersonl2007l) . The ab- 
solute astrometry is tied to the Guide Star Catalogue, version 
2.3.2 and has an uncertainty of 0.3". 

We use apertures of 3-pixel radius to measure fluxes and 
report the 1775 and 2^50 photometry for each epoch of each 
SN in the electronic version of this paper An example is 
shown in Table |2] The background noise is found empiri- 
cally by randomly placing apertures within regions that are 
free of objects and then measuring the dispersion in the inte- 
grated counts within these apertures. We have compared the 

' |http : //gsss ■ stsci ■ edu/Catalogs/GSC/GSC2/gsc23/gsc23. 



signal-to-noise ratios obtained for aperture and PSF photom- 
etry and found that the difference between the two is small 
if we use such small apertures. We use aperture photometry 
here, as it allows for a more robust correction of CTE and red 
halo scattering (described below). 

The photometry is corrected for variable CTE and for flux 
that is outside the aperture (referred to here as the aperture 
correction). CTE depends on the position of the source, the 
level of the background sky and the flux of the object that is 
being correcte d. It al so degrades with time. We follow the 
formulation of lRiess & Mack (2004) and apply the correction 
factor that corresponds to the 3-pixel radius case and the dates 
of our observations. On average, we applied a 4.3% correc- 
tion to the SN la flux. When calculating CTE corrections, we 
include the local background that had been previously sub- 
tracted. 

We also apply a color-dependent aperture correction on the 
ACS zgso-band measurements. In ACS zsso-band data, long- 
wavelength photons scatter off the back-side of the CCD, 
causing a degradation in the PSF. The effect is commonly 
known as red-halo scattering (ISirianni et al.lll998l) . The Tiny- 
Tim^ PSF does not account for this effect, hence it does not 
reproduce the observed zgso-band PSF. We studied how the 
PSF changes with wavelength using standard stars taken with 
a series of narrow -band filter observations from HST calibra- 
tion programs PID9020 and PID10720. We measured the red- 
halo scattering correction factor as a function of aperture ra- 
dius and wavelength. For a given aperture, we then can treat 
the correction as a modification of the F850LP throughput and 
zeropoint. We discuss the details of this procedure in the Ap- 
pendix. 

We have used the updated STScI ACS Vega zeropoints for 
light curve fitting. The latest zeropoints are from the STScI 
web site"* which were posted on May 19th, 2009. The ACS ze- 
ropoints changed on July 4th, 2006 due to a change in the de- 
tector temperature. The STScI definition assigns Vega (a Lyr) 
a magnitude of 0.0 in every filter However, in the Landolt 
system, which the SN la photometry in the literature refers 
to, Vega is not zero. To be consistent with literature SNe la, 
we introduce this non-zero Vega magnitude correction to our 
zeropoints. F or both the ^77.^ and z ssn-bands, the correction 
is 0.024 mag (iFukugita et al.lll996|). We therefore a djust the 
STScI zeropoints by this amount (lAstier et aDl2006l) . and use 
Vega zero points of ^775 = 25.291(STScI)+0.024 = 25.315 
and Z850 = 24.347(STScI)+0.024 = 24.371 for data taken 
before July 2006 (corresponding to a detector temperature of 
-77C), and ^775 = 25.277(STScI)+0.024 = 25.301 and zgso 
= 24.323(STScI)+0.024 = 24.347 for data taken after July 
2006 (corresponding to a detector temperature of — 81C). 

3.2. NICMOS Processing and Photometry 

All NICMOS science frames were processed with the lat- 
est CALNICA pipeline ( version 4.4. 1 iDahlen etaLll2008l) and 
then corrected for three well-known anomalies: the offset be- 
tween amplifiers, which affects all NICMOS exposures and 
is removed using the STSDAS PyRAF task PEDSKY; per- 
sistence after passage of the telescope through the South At- 
lantic Anomaly (SAA); and persistence after exposing the de- 
tectors to the limb of the Earth. Nine exposures are affected 
by the SAA, which leaves persistent signals from SAA cos- 
mic rays. We applied the STSDAS PyRAF task SAACLEAN 

^'http : / /www .stsc i ■ edu/ softwar e/tinytim/tinytim . h tml| 
_rel'^ila1steB.n/!);t.ffi;ffiwht-tii3ci ■ edu/hst/acs/ analysis/ zeropoints| 
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Table 2 

Photometry Data 



SN name 


Instrament 


Filter 


MJD 


Flux" 


Flux Error" 


Vega Zeropoint 


Exptime 


Nexp'' 


Raw Flux'^ 


Raw Flux Error*- 










(counts/s) 


(couiits/s) 




(s) 




(counts/s) 


(counts/s) 


SCP05D0 


ACS 


F850LP 


53564.098 


-0.0283 


0.0547 


24.371 


2000 


4 


-0.0188 


0.0365 


SCP05D0 


ACS 


F85()LP 


53589.117 


0.7733 


0.0573 


24.371 


2000 


4 


0.5153 


0.0381 


SCP05D0 


NICMOS 


Fl low 


53604.074 


0.5092 


0.0230 


23.029 


2560 


2 






SCP05D0 


ACS 


F850LP 


53610.836 


1.5084 


0.0649 


24.371 


2000 


4 


1.0040 


0.0432 


SCP05D0 


ACS 


F850LP 


53633.184 


0.7290 


0.0704 


24.371 


1500 


4 


0.4764 


0.0460 


SCP05D0 


ACS 


F775W 


53633.215 


0.7989 


0.2229 


25.315 


375 


1 






SCP05D0 


ACS 


F85()LP 


53654.434 


0.2388 


0.0694 


24.371 


1500 


4 


0.1520 


0.0442 


SCP05D0 


ACS 


F775W 


53654.469 


0.3908 


0.1725 


25.315 


375 


1 






SCP05D0 


ACS 


F850LP 


53679.258 


0.1792 


0.0708 


24.371 


1500 


4 


0.1137 


0.0450 


SCP05D0 


ACS 


F775W 


53679.273 


-0.1432 


0.1823 


25.315 


375 


1 






SCP05D0 


ACS 


F850LP 


53704.266 


0.0279 


0.0715 


24.371 


1500 


4 


0.0174 


0.0446 


SCP05D0 


ACS 


F775W 


53704.305 


0.0624 


0.2072 


25.315 


375 


1 






SCP05D0 


ACS 


F850LP 


53965.207 


0.0356 


0.0738 


24.347 


1360 


4 


0.0216 


0.0448 


SCP05D0 


ACS 


F775W 


53965.246 


-0.0034 


0.1981 


25.301 


515 


1 







Note. — The complete set of SNe la photometry data is available in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown 
here for guidance regarding its form and content. 

" Flux, corrected for CTE and color-dependent aperture correction (for ACS) and the count rate non-linearity (for NICMOS). For the ACS F850LP filter. 
zg50, this is the result of the iterative Method I in Appendix IaI 
The number of exposures. 

CTE corrected but not aperture-corrected flux for the ACS F85 OLP data. This flux column is used as an input for the modified filter Method II in Appendix 
IaI Note we use the modified filter response curve and shifted zeropoint as described in Appendix IaI 
The zeropoint has changed slightly after July 4th, 2006 (MJD=53920) due to the change in detector temperature. 

(iBarker et al.ll2007"l) to remove SAA persistence effects from 
the images. When a NICMOS observation is immediately 
preceded by an ACS data dump, the process could delay the 
NICMOS placement of the filter blank, subjecting the detec- 
tors to the bright limb of the Earth, whi ch imprints a per- 
sisten t pattern on subsequent exposures (iRiess & Bergeron] 
l2008h . Four exposures were affected in this way and were 
corrected using the STSDAS software NIC_REM_PERSIST. 
At this point, the mode of the flux distribution in each image 
is measured and recorded. These values are used as the sky 
levels for the count-rate non-linearity correction. 

Even after correcting NICMOS data for these well- 
studied anomalies, significant large-scale background non- 
uniformities remain. We developed methods to extract and 
remove the background structures; these are detailed in 
iHsiao et al] ( 1201 Ih Briefly, the models for the background 
structures are studied and characterized using approximately 
600 NICMOS exposures observed through the FUOW filter 
and processed with the procedures described above. Princi- 
pal component analysis applied on these images revealed that 
the intensity of the residual comer amplifier glow depends on 
the exposure sequence. The amount of residual glow decays 
exponentially and resets every orbit. With exposure times on 
the order of 1000 seconds, the exposures can be separated 
into two glow groups, each with approximately constant in- 
tensity. This makes it possible to extract the residual glow al- 
gebraically. The structured background is modeled as a com- 
bination of a constant component and a component that scales 
with the sky level and exposure time. The models are derived 
from the algebraic manipulation of stacked images for each 
glow group. The resulting constant component of the model is 
dominated by residual amplifier glow at the corners and resid- 
ual persistence structure at the center The model component 
which scales with sky level and exposure time displays a curi- 
ous fringe pattern whose origin is unknown. The model com- 
ponents are fit to individual exposures via scale parameters 
to create the customized background models to be subtracted 
from the individual exposures. In a final step, the bias offsets 



apparent in the middle column an d middle row are r emoved. 
Additional details can be found in lHsiao et alj ( 1201 Ih . 

3.2.1. NICMOS Count-rate non-linearity 

The NICMOS data are critically important for measuring 
the color of z > 1 SNe la. Any uncertainty in the NICMOS 
calibration severely limits the usefulness of SNe la observed 
with NICMOS. In particular, the NIC2 detector exhibits a 
count-rate dependent non-linearity (Bohlin et al. 2005), the 
severity of which is a function of wavelength. 

This non-linearity previously has only been studied at coun t 
rates three orders of magnitude higher (Ide Jong et al.l 120061 ) 
than the count rate of a typical SN la at 2; = 1, meaning 
that correcting the flux of SN la at z = 1 requires signifi- 
cant extrapolation and has a level of uncertainty that is diffi- 
cult to quantify. For example, comparisons between ground- 
based near-IR data and a different NICMOS camera (NIC3) 
showed that little or no co rrection is required for that camera 
dMobasher & Riessll2005h. It is difficult to r eco ncile this find- 
ing with the findings of lBohlin et al.l (l2006f) and lde Jong et alJ 
(2006)- A simple test at the flux levels relevant for the su- 
pemovae in this paper shows a difference of 12% between 
NIC2 and NIC3 when the non-linearity corrections are made, 
revealing significant problems with these extrapolations. 

For the NIC2/F110W filter, the degree of count-rate non- 
linearity is ~ 0.06 mag per factor of 1 change in count rate 
(iBohUn et alJ2006l:lde Jong et al.l2006h . The count rates from 
stars that are used to determine the NICMOS zero points are 
five orders of magnitude higher than the count rate of a typical 
SN la at 2 = 1. This corresponds to a ^ 0.3 mag correction 
for the NIC2/F110W filter! 

Since this is so important to the cosmological results, we 
have developed a method to address this count-rate non- 
linearity calibration directly (Ripoche et al. 2011). We an- 
alyze ACS, NICMOS, and ground-based near-IR observa- 
tions of early-type galaxies from clusters RCS J022 1.6— 0347 
(z = 1.02), RDCS J1252.9-2927 (z = 1.24), and XMMU 
J2235.3— 2557 (z = 1.39). The space and ground-based 
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data are used to constrain the spectral energy distributions 
(SED) of these galaxies, which are then numerically inte- 
grated through the Fl lOW filter transmission curve and com- 
pared to the counts measured with NICMOS. The principle 
advantage of the technique is that the count rate from early 
type galaxies at this redshift is similar to that measured for 
SNe la, i.e. about 0.03 counts/second/pixel (the contribution 
from amplifier glow is comparable). We applied this tech- 
nique using three galaxy clusters that have deep ground based 
near-IR imaging data from the VLT and deep images with the 
ACS and NICMOS camera. All three clusters are at different 
redshifts and produced consistent results. At the low count 
rates that are applic able to h igh-re dshift SNe la, we find tha t 
the prescription of iBohlinI (120071) and Ide Jong et all (12006b 
over-predicts the zeropoint correction for the NIC2 camera 
with the FllOW filter by 0.065 mag. We therefore use our 
zeropoint of 23.029 (Vega magnitude ) or 23. 757 (AB mag- 
nitude s). Additional details can be found in iRipoche et alj 
(l20Tll) . 

At high count rates, the count-rate non-linearity size 
has a strong depe ndence with wavelength across filters 
dde Jong et alj|2006b . being considerably stronger in bluer fil- 
ters. The SED of an early-type galaxy at z ^ 1.2, is a good 
match to a SN la about 20 rest-frame days after maximum, but 
is redder than a supernova SED at maximum (though this is 
compensated somewhat by the fact that the background level 
is about 1/3 of the source flux and is blue in the Fl lOW band- 
pass). The size of the count-rate non-linearity correction will 
thus also depend weakly on the phase and redshift, varying 
from 0.02 magnitudes at maximum to no additional correction 
20 rest-frame days after maximum. Since the wavelength- 
dependence of the non-linearity may not be even this strong 
at low count rates, we apply half the correction applicable at 
each phase, and add (in quadrature) an additional 0.01 mag- 
nitudes to the Fl low zeropoint error to account for this un- 
certainty. When added to the 0.006 mag statistica l error, and 
0.021 mag systematic error (IRipoche et al.ll20TTl) . this gives 
a total uncertainty on the zeropoint of 0.024 magnitudes. For 
the GOODS supernovae with N ICMOS o bservations, we start 
with the original flux given bv lRiess eTa l. (2007) (after con- 
verting the magnitude measurements to fluxes using the given 
zeropoint of 22.92), but increase the flux by 0.01 magnitudes, 
representing half the correction for the (possible) wavelength- 
dependence of the count-rate non-linearity. 

3.2.2. Galaxy Models 

Aft er the postprocessing described above (iHsiao et al.l 
we measure fluxes from the eight SNe la with NIC- 
MOS observations by performing PSF photometry on the im- 
ages. In all cases, the SNe la are not separated enough from 
their hosts to allow us to fit for the supernova flux alone; rather 
we fit a model of the host galaxy as well. By performing 
PSF photometry using analytic galaxy models, we avoid re- 
sampling the images (the better PSF sampling for the ACS 
data negates this advantage of PSF photometry), and extract 
the maximum possible signal-to-noise from our observations. 
We fit an analytic model of the host galaxy even when we 
have reference images, as this gives higher signal-to-noise, 
and nearly uncorrected photometry between epochs'. 

' Had we subtracted the flux in the reference images at the location of 
the supernova, the errors from this flux would have to be propagated as a 
covariance for all the other epochs. The errors on the galaxy model at the 
position of the supernova are typically much smaller 



Model PSFs for the supernovae are obtained with the 
TinyTim software using supernova SED templates from 
iHsiao et al] ( l2007h redshifted to the supernova redshift and 
warped as a function of wavelength to match the photometry. 
After the lightcurve fitting is complete, new PSFs are gener- 
ated from the SEDs based on this photometry and the process 
is repeated. Model PSFs for the galaxies are obtained with 
TinyTim by appropriately redshifting a galaxy spectrum from 
iBruzual & Charlog (120031) with an age of 2.5 Gyr and a solar 
metallicity; the exact shape of the galaxy spectrum does not 
greatly affect the results. The PSFs used are 3" in diameter, 
comparable to the patch fit in each NICMOS image. 

Although there is virtually no information at scales smaller 
than about half a pixel, all PSFs are seven times oversam- 
pled. This oversampling is necessary because the PSF is made 
slightly wider by the convolution with the subsampled pixels, 
increasing the flux of the derived photometry. In order for 
this effect to be negligible, seven times oversampling must be 
used. Finally, a correction is made to match the photome- 
try from the 3" TinyTim PSFs to the 30" TinyTim PSFs used 
in Ripoche et al. (2011). These differently-sized PSFs show 
different structure far in the wings, but the flux in the core 
changes by 3.5%, with negligible variation. 

We generally model the host galaxies as ellipsoids, with 
radial profiles given by second degree polynomial splines. 
These splines have ten nodes, with spacing that asymptot- 
ically approaches an exponential away from the core. The 
higher node density near the core provides more freedom to 
model the host where the flux changes quickly with position. 
In the few pixels closest to the core, where the spline changes 
rapidly, we numerically integrate over each subpixel befor e 
convolving with the PSF. On the basis of our tests (see i l3.2.3l ), 
the hosts of some supernovae were modeled with modifica- 
tions to this basic scheme, as discussed in the following sec- 
tion. 

3.2.3. Photometry Testing 

Three ingredients all have to be correct in order to achieve 
photometry with low bias and variance: the PSF model, the 
galaxy model, and the supernova centroid. Deriving a PSF 
from a field star (details in the SN SCP06C0 discussion be- 
low) and comparing against TinyTim gives photometry con- 
sistent to a few mmags, so we do not believe this is a major 
contribution to our errors. 

Testing the host galaxy model and supernova centroiding is 
more involved. For each observation, we subtract the best-fit 
supernova light, and place simulated supernovae (at the same 
flux level) in the images. The only place one cannot do this 
test is at the location of the actual supernova, as putting a sim- 
ulated PSF in this location yields a measurement that will be 
highly correlated with the measurement of the supernova. We 
therefore do not place any simulated supernovae closer than 
two pixels to the best-fit location of the supernova. By exam- 
ining the bias and variance of the extracted fluxes from a large 
number of simulations (~ 100), we can choose the galaxy 
model which gives the most precise and accurate^ fluxes for 
each particular supernova. We emphasize that the results of 
these simulations were the only metric used in choosing the 
detailed model. In particular, there was no feedback from the 
shape of the lightcurve or the Hubble diagram since these 
would have undercut the principles of "blind" analysis we 

^ We found that precision and accuracy coiTelated in our simulations. 
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tried to maintain (see The same basic galaxy model (dis- 
cussed above) was used for the NICMOS photometry of each 
supernova, with the following exceptions. 

• SN SCP06C0: As mentioned in ^TT] there is a small 
galaxy about 0.6" from the likely host of SN SCP06C0, 
and just 0.2" from SN SCP06C0 itself. We note that 
the surface brightness of the small galaxy is one fourth 
of that host at the location of the supernova. The host 
also has some azimuthal asymmetry visible, indicating 
a possible merger The cluster XMMl 229+01 was also 
observed as part of a program to cross-calibrate NIC- 
MOS (Ripoche et al. in preparation) and deep, well- 
dithered images were obtained in the WFC3 Fl lOW fil- 
ter, allowing a more-flexible background model to sub- 
tract both galaxies. We modeled the galaxies with a 
2D second-order spline, with nodes placed in a grid ev- 
ery 0.076" (the natural pixel scale of NICMOS). The 
WFC3 FllOW PSF was modeled as a combination of 
the elliptical galaxy model and a 2D spline (with a spac- 
ing of 0.1") using dithered images of a field star (This 
is the same empirical PSF model used for testing Tiny- 
Tim for NICMOS, although there the 2D spline nodes 
are spaced at the natural pixel scale of NICMOS.) Our 
testing indicates that this method achieves the same 
signal-to-noise ratio as the other supernovae that have 
simpler galaxy subtractions. 

• SN SCP06A4: We found a small amount of azimuthal 
asymmetry in the host. Adding a second-order 2D 
spline to the galaxy model, with a node spacing of 
0.38" (five times the natural pixel scale of NICMOS 2) 
successfully modeled this asymmetry, without adding 
additional measurement uncertainty to the supernova 
flux. 

• SN SCP05D6: The host galaxy requires two elliptical 
components to be fitted well. These components are 
forced to have the same centroid, but are allowed dif- 
ferent orientations, ellipticities, and radial profiles. One 
component forms a bulge, while the other one forms a 
disk. In one epoch contaminated by the SAA, aper- 
ture photometry with a one-pixel radius aperture on the 
galaxy-model-subtracted images gave better signal-to- 
noise than PSF photometry, so we used this instead. 

• SN SCP06U4: This supernova was on the core of a 
galaxy that appears to be merging with another galaxy. 
Similarly to SN SCP05D6, a second efliptical com- 
ponent was needed to model the host, (in this case, 
a third, detached component was used to model the 
fainter companion). Our simulated supernovae revealed 
that, rather than using one host galaxy model to extract 
photometry, even more precise results were obtained 
averaging photometry results derived using the ellipti- 
cal model and the 2D spline model (discussed above for 
SN SCP06C0). Using this procedure results in a change 
in flux well inside the error bar 

• SN SCP06H5: The one NICMOS observation of this 
supernova was our most challenging extraction. The 
observation of the supernova was 11 rest-frame days 
after maximum, and it is only ^ 0.1" from the core. 
As with SN SCP05D6, the host galaxy requires two el- 
liptical components to be fitted well. (Comparing to the 



2D spline model discussed above, we obtain photom- 
etry that is the same to within a small fraction of the 
error bar) 

The signal-to-noise ratio of this measurement is low, 
likely implying some amount of bias due to centroiding 
error However, this is the only measurement with a 
signal-to-noise this low, so no correlation is introduced 
with any other measurement. 

3.3. Keck AO Photometry 

The photometry of 2 > 1.2 SNe ha s been almost ex- 
clusively measured from HST images (|Knop et alj 120031 : 
iRiess et alJl2004l2007HAmanullah et al.ll20IOh . At these red- 
shifts, the rest frame _B-band is redshifted beyond 9000 A, 
and falls in the near-IR (NIR). NIR observations are typically 
much easier from space, not only because of the higher spatial 
resolution of HST compared to ground-based seeing-limited 
systems, but also because of the lower NIR sky noise in space. 
However, we show here that adaptive optics on large ground- 
based telescopes can overcome these limitations and allow 
high-z SNe to be studied from the ground. 

We observed the z=1.315 SN SCP05D6 with the Keck 
Laser Guide Adaptive Optics (LGS AO) system. These ob- 
servations were made in the i/-band (1.6 /im), which corre- 
sponds to the rest-frame i?-band. The diffraction-limited res- 
olution of the Keck AO images was ^ 0.05", or a factor of 
three better than the spatial resolution of HST at these wave- 
lengths. The high spatial resolution meant a much better sep- 
aration of the SN from the galaxy core compared to HST. It 
also allows greater contrast between the SN and the sky back- 
ground. We obtain a photometric precision of 0.14 mag at 
H^24 mag in a one hour exposure with Keck AO, showing 
the potentia l of AO in SN la cosmology. 

Melbour ne et all (^2007^ reports the details of the Keck AO 
photometry, here we briefly summarize the observations. The 
Keck LGS AO observ ing runs for the Center f or Adaptive Op- 
tics Treasury Survey ([Melbourne et al.l [2005b coincided with 
the HST Cluster SN Survey program, and we successfully ob- 
served SN SCP05D6 at three epochs, before, near and after 
the lightcurve maximum. 

A 14 mag star, 25" away was chosen to provide AO tip-tilt 
correction. The ~ 11 mag sodium LGS was pointed at the 
galaxy to provide higher-order AO corrections. The obser- 
vations were sampled with a 0.01" pixel-scale, allowing the 
diffraction-limited Keck PSF to be fully resolved. Individual 
exposures of 60s were taken with five dithered pattern posi- 
tions. The sequence was repeated until sufficient depth was 
reached. Total exposure times varied from 30 min to 1 hour 
per epoch. 

We were also fortunate to have a 17.9 magnitude natural 
PSF star only 4" away from SN SCP05D6, so the PSF near 
the location of the SN was well-determined. From the star, 
we measured a FWHM of 0.055" while we had a mean H- 
band seeing of 0.4". Using the observed PSF , the host galaxy 
was modeled by GALFIT jPenget all 120021) and subtracted 
from the image, providing a clean measurement of the SN 
diffraction-limited core. Relative photometry with respect to 
the nearby PSF star was performed at each epoch, and cal- 
ibrated by the UKIRT standard star FS6 (the photometry is 
reported in Table |2). The photometric uncertainty was esti- 
mated by simulations of model PSFs embedded into the AO 
image at the same galacto-centric radius as the actual SN. 

We fit the SN-Hghtcurve with the photometric data from 
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HST/ACS F775W, F850LP, HST/NICMOS FllOW, and this 
Keck/ AO observation. We found that the Keck AO data was 
consistent with the HST observations (Figure |2]i. Although 
the uncertainty in the AO measurement was larger than that of 
HST, including it reduces the color uncertainty by more than 
10%, and reduces the sensitivity of the fit to the NICMOS 
zeropoint uncertainty by more than a third. 

4. AUGMENTING THE UNION2 SUPERNOVA 
COMPILATION: UN10N2.1 

SNe la are an excellent probe of dark energy, as they mea- 
sure the magnitude-redshift relation with very good precision 
over a wide range of redshifts, from z = 0uptO2~1.5 and 
possibly beyond. While some individual sets of SNe la are 
now, by themselves, large eno ugh to provide constraints on 
some cosmological parameters (iGuy et al.l2010l:lKessler et aP 
120091) . they do not yet constrain the properties of dark energy 
as well as analyses that combine individual data-sets to create 
a compil ation of SNe la t hat co vers a broader range of red- 
shifts. In iKowalski et alj ( 120081) . we developed a systematic 
methodology for combining the many available datasets into 
one compilation, called the "Union" compilation. 

There are many positive features behind the philosophy 
adopted by the Union analysis. It includes all SN la data- 
sets on an equal footing, with the same lightcurve fitting, 
cuts, and outlier rejection. Estimates of the systematic er- 
ror are entered into a covariance matrix, which can be used 
for fitting any cosmological model. Choices about how to 
do the analysis and what cuts to apply are done with the 
cosmological results hidden. This type of "blind" analy- 
sis mitigates biases that arise fr om inadvertently scrutin izing 
some data more than others. In lAmanullah et all (120101) . we 
adopted this strategy to create the Union2 compilation. This 
paper also revised and improved the Union analysis in sev- 
eral significant ways. Firstly, it augmented the Union sam- 
ple with new SN la data-sets from the literatu re, including 
102 lo w-redshft SNe la from the CfA3 survey (iHicken et al.l 
I2009al) . 12 9 intermediate-redshif t SNe la from the SDSS 
SN survey (iHoltzman et al.ll2008h . five intermediate-redshif t 
SNe la discovered from La Palma (lAmanuUah et al.l 120081) . 
and six new high-redshift SNe la. The paper revised the 
analysis by replacing the SALT lightcurve fitter with SALT2 
dGuy et al. 2005, 2007), and handled many systematic errors 
on a supernova-by-supernova basis in a covariance matrix. 

In the current paper, we use the analysis procedure that 
was used for the Union2 compilation with only one signifi- 
cant change: a correction for the host-mass SN la-luminosity 
relation, described below. The HST calibration and the asso- 
ciated errors have also been updated, as described in Section 
14.4.11 We refer to this new compilation as "Union2.1." 

4. 1 . Host Mass Correction to SN la Luminosities 

There is evidence that SN la luminosity correlates with 
the mass of the host galaxy, even after the cor rections for 
color and light curve width have been applied (iKellv et al.l 
120 lot ISullivan et aLll2010t iLampeitl et alllioiOl) . Since low- 
redshift SNe la are predominantly from surveys that target 
catalogued galaxies, the host galaxies of SNe la in these sur- 
veys are, on average, more massive than the host galaxies of 
distant SNe la from untargeted surveys. SNe la from low- 
redshift samples therefore have brighter absolute magnitudes. 
Left uncorrected, the correlation biases cosmological results 
dSuUivan et al..2010.) . 



ISulUvan et al.l ( 120101) find that the correlation can be cor- 
rected by fitting a step in absolute magnitude at tjj threshold _ 
lO^^mQ. There are two complications with making this cor- 
rection: most of the SNe in the Union2 compilation do not 
have host mass data available in the literature, and SN la hosts 
with masses close to the cutoff may scatter across, decreas- 
ing the fitted size of the step. To address these problems, we 
adopt a probabilistic approach to determining the proper host 
mass correction to apply to each supernova, correcting each 
supernova by the probability that it belongs in the low-host- 
mass category. (The low-host-mass category was chosen be- 
cause most of the low-redshift supernovae are from high-mass 
galaxies, so correcting the low-host-mass supernovae mini- 
mizes the correlation between Mb and the correction coeffi- 
cient.) 

Suppose we have a mass measurement m"*^^ and we would 
like to estimate the probability that the true mass m}^^'^ is less 
than the mass threshold. We begin by noting that 



p(TOf^mr'=) = 



(1) 



We can then integrate this probability over all true host masses 
less than the threshold: 



/'(mf"'= < '■'="'^°''^|mf = (2) 

^threshold 

P(mf"|mf"'=)P(mf"°) 



up to a normalization constant found by requiring the inte- 
gral to be unity when integrating over all possible true masses. 
P(m^'"°) is estimated from t he observed distribut ion for each 
type of survey. The SNLS (ISullivan et al.1 120101) and SDSS 
(|Lampeitl et al. 2010) host masses were assumed to be rep- 
re sentative of untarg eted surveys, while the mass distribution 
m IKellv et alJ (iloTol) was assumed typical of nearby targeted 
surveys. As these distributions are approximately log-normal, 
we use this model for P(m*'"°) using the mean and RMS 
from the log of the host masses from these surveys (with the 
average measurement errors subtracted in quadrature), giv- 
ing logio -P("^l™°) = A/'(m = 9.88, cr2 ^ 0.922) 
targeted surveys and log^Q P(m*™'=) = 7V(10.75, 0.ee^) for 
targeted surveys. When host mass measurements are avail- 
able, P(mf "|m*™° ) is also modeled as a log-normal; when 
no measurement is available, a flat distribution is used. 

For a supernova from an untargeted survey with no host 
mass measurement (including supernovae presented in this 
paper which are not in a cluster), P(m^''"° < threshold^ 
is the integral of P(7n^'"°) up to the threshold mass: 0.55. 
Similarly, nearby supernovae from targeted surveys with- 
out host galaxy mass measurements are given a P(to^''"° < 
^threshold) (ygj.y similar numbers of 0.50 and 0.09 

are derived from the observed distribution, without using the 
log-normal approximation.) We must make the connection for 
supernovae in clusters, as these are from a targeted survey. We 
take advantage of the simpler SEDs of early-type galaxies to 
precisely measure these masses ' . 

The best-fit mass-correction coefficient, 5, is much smaller 
in magnitude (—0.03) than that found in other studies (« 

' C-OOl and F-012 are in clusters, but are not hosted by early-type hosts. 
We use the untargeted value for their host-mass-luminosity relation correc- 
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—0.08). This may be due to the small value for 5 from the 
first-year SNLS data, as shown in Table |6] We in clud e the 
difference in Ss as a systematic, as discussed in ^4.5\ For 
this analysis, we assumed the host-mass correction does not 
evolve with redshift. 

4.2. Light-Curve Fitting 

Following A manullah et all (120101) . we use SALT2 
(iGuy et al. 2007^ to fit supernova lightcurves. The SALT2 
model fits three parameters to each SNe: an overall normal- 
ization, xq, to the time dependent spectral energy distribu- 
tion (SED) of a SN la, the deviation, xi, from the average 
lightcurve shape, and the deviation, c, from the mean SN la 
B — V color The three parameters, xi, c, and integrated B- 
band flux of the model SALT2 SED at maximum light, mg^^, 
are then combined with the host mass to form the distance 
modulus 



(3) 

where Mb is the absolute B-band magnitude of a SN la with 



0, c = and F(to^™'= < m 



threshold 



) = 0. The param- 



Xi 

eters a, {3, S and Mb are nuisance parameters that are fitted 
simultaneously with the cosmological parameters. The SN la 
photometry data and SALT2 light curve fits are shown in Fig- 
ure |2] The fitted S ALT2 parameters are listed in Table |3] as 
well as the host galaxy host stellar mass and lensing magnifi- 
cation factor. 

4.3. Union!. 1 

To the Union2 SN la compilation ( lAmanullah et aljr2010l) . 
we add 16 SNe la from this paper that were classified as ei- 
ther secure or probable, including six SNe la hosted by high-z 
cluster elliptical galaxies. The four SNe la that were classi- 
fied as possible are not u sed. We also add 18 S Ne la from 
the low -redshift sample of IContreras et al.l (120101) . 9 of which 
were not in Union2 (the others had published data from CfA). 
As in Union2, for all SNe we require 

1. that the CMB-centric redshift is greater than 0.015; 

2. that there is at least one point between —15 and 6 rest- 
frame days from _B-band maximum light; 

3. that there are at least five valid data points; 

4. that the entire 68% confidence interval for xi lies be- 
tween —5 and +5; 

5. data from at least two bands with rest-frame central 
wavelength coverage between 2900 A and 7000 A; and 

6. at least one band redder than rest-frame U -band (4000 
A). This cut is new to this analysis, but only affects 
SN 2002fx, a GOODS supernova which is very poorly 
measured. 

In addition to these quality cuts, we removed any supernova 
spectroscopically classified as SN 1991bg-like. These SNe la 
are a distinct subclass which is not modeled well by SALT2. 
At high redshift, where spectroscopic sub-typing may not be 
possible, we screen for these supernovae photometrically by 
searching for any supernovae with red (c > 0.2) and narrow- 
width {xi < —3) lightcurves, but do not find any. When 



fit with SALT2, and color-corrected and shape-corrected (as 
though they were normal SNe la), spectroscopically identi- 
fied members of this class have an average absolute magnitude 
only 0.2 magnitudes fainter than normal SNe la; any contam- 
ination from the handful of supernovae near this cut will have 
only a small impact (and one well-accounted for b y our con- 
tamination systematic, see lAmanullah et aTl (120101) ). 

From the 16 SNe la that were classified as either secure 
or probable (see Table SN SCP06U4 and SN SCP06K18 
fail to pass these cuts. SN SCP06K18 lacks good enough 
light-curve coverage and SN SCP06U4 fails the xi cut'. This 
leaves 14 SNe la that are used to constrain the cosmology. 

4.4. Fitting the Cosmology 

Following lAmanullah et al.l (120101) . the best-fit cosmology 
is determined by minimizing 



Astat 



SNc 



cr, 



sample 



A detail ed discussion of t he ter ms in this equation can be 
found in lAmanullah et all (120101) . We only comment on the 
final term in the denominator, CTg^j^pie, which is computed by 
setting the reduced of each sam ple to unity. This term 
was referred to as "fgygj^j^ajij," in iKowalski etal] ( 120081) ; 
lAmanullah et al.l (120101) . We note that o'^^^^yc includes intrin- 
sic dispersion as well as sample-dependent effects. This term 
effectively further deweights samples with poorer-quality data 
that has sou rces of error which have not been accounted for 
As noted in lAmanullah et al.l ( 120101) . this may occasionally 
deweight an otherwise we l l-measured supernova. 

Followin g IConlev etal J ( 120061) . IKowalski etall ('2008) and 
lAmanullah et al .1 (120101) . we hide our cosmology results un- 
til the full analysis approach is settled. As in previous Union 
analysis, we carry out an iterative minimization with out- 
lier rejection. Each sample is fit for a flat ACDM cosmology 
independently of the other samples (but with a, j3, and 5 set 
to their global values). An Mb is chosen for each sample 
by minimizing the absolute variance-weighted sum of devia- 
tions, minimizing the effects of outliers. We then reject any 
supernova more than 3(t from this fit. All of the SNe la in 
our new sample pass the outlier rejection. As each sample is 
fit independently with its own Hubble line, systematic errors 
and the choice of cosmological model are not relevant in this 
selection. 

4.4.1. Diagnostics 

A diagnostic plot, which is used to study possible incon- 
sistencies between SN la samples, is shown in Figure [3] The 
median of (isampic can be used as a measure of the intrinsic 
dispersion associated with all SNe la. The intrinsic dispersion 
is a reflection of how well our empirical models correct for 
the observed dispersion in supernova luminosities. The me- 
dian (Jsampic for this paper is 0.15 mag and is indicated with 
the leftmost dashed vertical line in the left panel. 

The variance weighted RMS about the best-fit cosmology 
gives an indication of the quality of the photometry. A sample 
with more accurate photometry will have a smaller RMS. For 

' Using an updated version (2-18-17) of SALT2 (or using SALTl), 
SN SCP06U4 would pass this cut, so this supernova may be included in future 
analyses. 
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Observed Date : MJD - 53963.637 Observed Date ; MJD - 53718.263 Observed Date ; MJD - 53735.334 
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Observed Date : MJD - 53848.254 Observed Date ; MJD - 53675.59 Observed Date ; MJD - 53840.789 
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Figure 2. 15 SNe la light curve fits by SALT2. Flux is nomialized to the 2g5o-band zeropoint magnitude. ACS 4775, ACS zgso and NICMOS Fl lOW data is 
color coded in blue, green and red respectively. Note that SCP 05D06 (z=1.314) has i/-band data from Keck AO system (orange) iMelboume et aL,2Q07ft and 
that this data is consistent with the HST/ACS and HST/NICMOS light curve data. 
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Table 3 

SALT2 Lightcurve Fit Results 



SN name 


z 


MJDBmax 






c 




Galaxy Mass" 
(10"Mq) 


Lens 
Factor'' 


SCP06A4 


1.192 


53912.7 ± 1.5 


25.497 ± 0.048 


-1.45 ±0.68 


0.065 ±0 


084 


0.44 




SCP06C0 


1.092 


53735.4 ± 1.0 


25.636 ± 0.066 


-2.66 ±0.65 


0.257 ±0 


083 


1.97 


1 030+°°°'^ 

l.UJU_Q ggg 


SCP06C1 


0.980 


53759.0 ±0.7 


24.613 ± 0.028 


-0.35 ±0.33 


0.014 ±0 


053 




SCP06F12 


1.110 


53718.4 ±2.3 


25.253 ±0.068 


-2.09 ± 1.29 


-0.133 ±0 


142 






SCP06G4 


1.350 


53860.9 ± 1.4 


25.424 ±0.052 


0.15 ±0.64 


-0.029 ±0 


052 


1.72 


i.U10_g 


SCP06H3 


0.850 


53848.2 ± 0.6 


24.345 ± 0.038 


0.58 ±0.31 


0.089 ± 


067 




SCP06H5 


1.231 


53860.2 ± 1.5 


25.389 ±0.111 


-3.12 ± 1.10 


-0.103 ±0 


187 


3.66 




SCP06K0 


1.415 


53751.3 ±2.8 


25.811 ± 0.087 


0.30 ±0.97 


0.147 ±0 


081 


2.30 




SCP06N33 


1.188 


53962.6 ±4.3 


25.407 ±0.132 


-2.15 ± 1.32 


-0.038 ± 


175 




1 066+°'^" 

i.UDO_Q gj^4 


SCP05D0 


1.014 


53606.9 ±0.9 


25.201 ± 0.066 


-0.61 ±0.65 


0.061 ±0 


085 


0.40 


SCP05D6 


1.315 


53658.5 ± 1.3 


25.660 ± 0.046 


-1.26 ±0.56 


-0.058 ±0 


061 


2.61 


1 021+°"12 


SCP05P9 


0.821 


53675.6 ±0.6 


24.367 ± 0.049 


0.25 ±0.50 


0.022 ±0 


075 




SCP06R12 




53966.6 ±3.5 


25.789 ±0.114 


-2.06 ± 1.50 


-0.158 ±0 


198 


0.23 




SCP06U4 


1.050 


53944.4 ± 1.1 


25.056 ± 0.063 


-4.62 ± 1.09 


-0.102 ±0 


096 


1.11 




SCP06Z5 


0.623 


53840.5 ± 3.0 


23.482 ± 0.144 


-0.76 ±0.88 


0.070 ± 


120 







" The details of host galaxy identifications, coordinates and its stellar mass measurements can be found in lMevers et al.l J201 II) . 
^ Gravitational lensing magnification factor (see i]2. ll for details). For cosmological analysis we must divide the corrected SNe fluxes 
by this factor to make use of these supemovae. 
SCP06U4 is not included in our current cosmological results, but will likely be included in future compilations (see [|4]for details). 
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Figure 3. Diagnostics plot for the individual data sets. From left to right: irreducible sample dispersion (filled circles) and variance-weighted RMS about the 
best-fit model (open circles); the average sample residual from the best-fit model (^measured ~ Mmodcl) excluding and including systematic errors; and the 
best-fit slope of the Hubble residual (in magnitudes) versus redshift — 9/Jrcsidual/9z. Note that the eiTors on the sample dispersion include only statistical 
errors and do not include possible systematic errors. The confidence intervals on the wei ghed R MS are obtained with Monte-Carlo simulations. The triangles in 
the sample residual plot show the effect of including the filter shifts discussed in Section (4.4. II 
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SNe la from our survey, the RMS is 0.19 ± 0.04, which is 
only slightly larger than that measured for the 1st year SN la 
sample from SNLS, and equal to the median of all samples 
(shown as the rightmost dashed line in Figure|3] left panel). 

The two middle panels show the tension between data-sets, 
the first with statistical errors only, and the second with sta- 
tistical and systematic errors (see ^4.5i . Most samples land 
within 1 (7 of the mean defined by all samples and about one 
third lie outside 1 a, as expected for a normal distribution. No 
sample exceeds 2 a. The right hand panel shows the slope of 
the residuals, which, for larger data sets, can be used to reveal 
Malmquist-like biases or calibration errors. 

The supernovae from our sample are 1 .5-a brighter than the 
average sample. While the source of the difference may cer- 
tainly be a simple statistical fluctuation, part of the difference 
might be attributable to errors in the filter responses of the 
ACS filters. (The difference is largely driven by the SNe la 
that have only ACS ^775 and 2350 data to constrain their light 
curves.) Based on photorn e tric ob servations of spectropho- 
tometric standards, iBohlinI (l2007h report possible blueward 
shifts of 94 A for the zgso filter and 57 A for the 1775 fil- 
ter (with smaller shifts in bluer filters). The red triangle in 
the sample residual panel shows the effect of applying these 
shifts. The shifts also affect the GOODS supernovae. The 
green trian gle shows the affect of applying the filter shifts to 
those data. IBohlinI (120071) notes that more data to confirm the 
filter shifts are needed, so we do not apply them in our pri- 
mary analysis. Instead, we include the uncertainty in the filter 
curves as a systematic error, as described in ^i5\ 

Part of the difference could also be due to the correction that 
we apply for the recently discovered correlation between host 
galaxy mass and the luminosity of SNe la after the lightcurve 
width and color corrections have been applied. Many of the 
hosts in our sample are massive early type galaxies. In this 
analysis, the correction we use is smaller than the correction 
that has been noted by others. We add this difference as a 
systematic error, as described in ^4.5\ 

Figure |4] shows the Hubble Diagram with SNe from the 
updated Union2 sample and the best-fit ACDM model. We 
add 14 SNe la from this current paper. (As discussed above, 
SN SCP06U4 is likely to be included in future analyses so it 
is included on the plot with a different symbol.) Ten (eleven 
with SN SCP06U4) are above a redshift of one, significantly 
increasing the number of well-measured supernovae above 
this redshift. 

4.5. Systematic errors 

In this paper, we follow the sy stematics analysis we pre- 
sented in lAmanuUah et all (120101) . Systematic errors that di- 
rectly affect supernova distance measurements (calibration, 
and galactic extinction, for example) are treated as nuisance 
parameters to be fit simultaneously with the cosmology. Min- 
imizing over these nuisance parameters gives additional terms 
to add to the distance modulus covariance matrix 



where the sum is over each of these distance systematic errors 
in the analysis. (Although the distance modulus depends on 
5 as well as a and /3, the derivatives with respect to the zero- 
points do not.) In this analysis, a and /3 have little interaction 
with cosmological parameters. When computing cosmologi- 
cal constraints, we therefore freeze the covariance matrix in 



order to avoid multiple matrix inversions . Only when the a 
and P may vary significantly from the global best-fit (Table 
|6]l, do we update a and /3. 

Systematic errors that affect sample composition or the 
color and shape correction coefficients cannot be parameter- 
ized supernova-by-supemova in this way. These are incorpo- 
rated by assigning each dataset its own constant covariance. 
This is an adequate treatment, as these systematic errors are 
subdominant. 

There are two systematic errors that were not included in 
lAmanuUah et al. (2010), but are included in this analysis for 
the first time: a systematic error on the host-mass correction 
coefficient, 5 (which might affect 6 at the ~ 0.05 level), and 
uncertainties in the effective wavelengths of the ACS ^775 and 
Z850 filters. 

In addition to updating the NICMOS Fl lOW zeropoint and 
uncertainty, as described in ij3.2.1l we revise the uncertainty 
assigned to the zeropoint for NICMOS F160W to account for 
the un certainty in the coun t-rate non-linearity at this wave- 
length dde Jong et al.ll2006l) . Table g] gives the assumed zero- 
point error for each filter. 

We note that the nearby supernovae from targeted searches 
are sensitive to 5 (relative to the untargeted searches) at the 
level of (0.55 - 0.13)AJ ?» 0.02 magnitudes, while the co- 
variance weighted mean of the cluster supernovae varies with 
5 as 0.24A(5 k, 0.01 magnitudes. We cannot propagate this 
systematic on a supernova-by-supemova basis, as this would 
be equivalent to fitting for 5, which we already do. Therefore, 
we include this error by adding a covariance of 0.02^ to the 
nearby, targeted supernova surveys, a covariance of 0.01^ to 
our new data-set, and 0.02 • 0.01 between these data-sets. 

Including uncertainties in filter effective wavelength is not 
as straightforward as including zeropoint uncertainties. Ef- 
fective wavelength is only the first-order method of describ- 
ing a filter. For a simple filter shift, as implemented here, 
(i/i(a, /3) / d\ will undergo significant variations as supernova 
spectral features shift in and out of the filter. These are likely 
to be worse than the actual effect of simply reweighting fil- 
ter throughput. Although in general these variations will get 
averaged out with different phases, redshifts, and additional 
filters, we have modeled a worst-case in accounting for this 
systematic (and even then it only affects the supernovae most 
dependent on zsso)- 

Table |5] shows the impact of each type of systematic error 
on our cosmological constraints, in combination with BAO, 
CMB, and Ho data (see IQ. For the purpose of constructing 
Table |5] we add, for each systematic error in the table, the 
contribution from just that systematic to the statistical-only 
covariance matrix. The confidence interval for constant w 
where the is within 1 of the minimum (the edges of this 
confidence interval are hereafter referred to with the notation 
= 1) is found iteratively; the plus and minus errors for 
constant w are averaged. The statistical-only constant w er- 
ror bar is subtracted in quadrature, leaving the effect of each 
systematic on constant w. We also quote the effect of each 
systematic error of the Ax^ = 5.99 confidence contour in the 
{wotWo) plane; as this is two-dimensional, we subtract the 
area (not in quadrature) of the statistical-only contour. 

Since the derived cosmology errors vary with the best-fit 
cosmology, after a given systematic error has been added, the 

' As demonstrated in the Union2 appendix, tliese matrix inversions can be 
simplified at the expense of more matrix multiplication; the run-time does not 
change much. 
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Cluster Search (SCP) 
Amanullah et al. (2010) (SCP) 
Riess et al. (2007) 
Tonry et al. (2003) 



Contreras et al. (2010) 
Hicken et al. (2009) 
Kowalski et al. (2008) (SCP) 

Riess et al. (1999) 

Hamuy et al. (1996) 



Miknaitis et al. (2007) 
Astier et al. (2006) 
Knop et al. (2003) (SCP) 
Amanullah et al. (2008) (SCP) 
Barris et al. (2004) 
Perlmutter et al. (1999) (SCP) 
Riess et al. (1998) + HZT 
Holtzman et al. (2009) 



0.2 



0.4 



0.6 0.8 

Redshift 



1.0 



1.2 



1.4 



Figure 4. Hubble diagram for the Union2. 1 compilation. The solid line represents the best-fit cosmology for a fiat ACDM Universe for supemovae alone. 
SN SCP06U4 falls outside the allowed x'l range and is excluded from the current analysis. When fit with a newer version of SALT2, this supemova passes the 
cut and would be included, so we plot it on the Hubble diagram, but with a red triangle symbol. 



Table 4 

Assumed instmmental uncertainties for SNe in this paper 



Source 


Band 


Uncertainty 


Reference 


HST 


WFPC2 


0.02 


Hever et al. (2004) 




ACS F850LP 


0.01 


Bohlin (2007) 




ACS F775W 


0.01 






ACS F606W 


0.01 






ACS F850LP 


94 A 


Bohlin (2007) 




ACS F775W 


57 A 




ACS F606W 


27 A 






NICMOS J 


0.024 


RiDOche et. al. (inDrep). Sectionl3.2. 11 




NICMOS H 


0.06 


de Jong et al. (2006) 


SNLS 


g. r, I 


0.01 


Astier et al. (2006) 




z 


0.03 




ESSENCE 


R,l 


0.014 


Wood^Vasev et al. (2007) 


SDSS 


u 


0.014 


Kessler et al. (2009) 




g, r, i 


0.009 




z 


0.010 




SCP: Amanullah et al. f2010) 


R,I 


0.03 


Amanullah et al. (2010) 


J 


0.02 


Other 


U-band 


0.04 


Hicken et al. (2009a) 




Other Band 


0.02 


Hicken et al. (2009a) 
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Table 5 

Effect on constant w eiTor bars and area of the 95% wq — Wa confidence contour (inverse DETF 
FoM) for each type of systematic error, when SN la constraints are combined with constraints 
from CMB, Hq, and BAO. 



Source 


Error on Constant w 


Inverse DETF FoM 


Vega 


0.033 


0.19 


All Instrument Calibration 


0.030 


0.18 


(ACS Zeropoints) 


0.003 


0.01 


(ACS Filter Shift) 


0.007 


0.04 


(NICMOS Zeropoints) 


0.007 


jO.Ol 


Malmquist Bias 


0.020 


0.07 


Color Correction 


0.020 


0.07 


Mass Correction 


0.016 


0.08 


Contamination 


0.016 


0.05 


Intergalactic Extinction 


0.013 


0.03 


Galactic Extinction Normalization 


0.010 


0.01 


Rest-Frame ?7-Band Calibration 


0.009 


iO.Ol 


Lightcurve Shape 


0.006 


jO.Ol 


Quadrature Sum of Errors/ Sum of Area (not used) 
Summed in Covariance Matrix 


0.061 
0.048 


068 
0.42 
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supernova magnitudes are shifted so that the best-fit cosmol- 
ogy including that systematic matches the best-fit with sta- 
tistical errors only. This magnitude adjustment (which is the 
same adjustment we use for blinding ourselves to the best-fit 
cosmology) consists of repeatedly computing the difference 
in distance modulus between the best-fit cosmology and fidu- 
cial value and adding it to the supernovae. 

As with the Union2 compilation, calibration systematics 
represent the largest contribution to the error on constant w. 
Here, we see that they are also the dominant systematics for 
(wq, Wo)- As noted by Amanullah et al. (2010), significantly 
smaller systematic errors are derived by adding each covari- 
ance in the covariance matrix, rather than adding the cosmo- 
logical impacts together This is due to the different red- 
shift dependence of each systematic error, as well as some 
self-ca libration that occurs as described in lAmanullah et al.l 
(l2OT0l) . 

Some potential systematic errors can be investigated by di- 
viding the whole dataset into subsets. Table |6] shows many 
of these divisions. All of the numbers are computed includ- 
ing supernova systematics; the cosmological constraints are 
computed including BAO and CMB data. In short, we do not 
see any evidence of unknown systematic errors, requiring the 
cosmological impact to be smaller than the current errors. 

The first subsets are subsets in redshift. These can be 
used to study possible evolution of correction coefficients for 
shape, color, and host mass. The redshift range 0.5 to 1 seems 
to show p and 5 smaller in magnitude, but the revised SNLS 
sample jGuv et al.ll2010l) which uses a newer version of the 
calibration and lightcurve fitting (as well as many more super- 
novae), shows no signs of this. As we have already budgeted 
these systematic uncertainties, these updates will be within 
our error bars. 

The next rows show the effect of changing 5 from to 
—0.08 (the size of the correction in Sullivan et al. (2010)). 
Because a large error on 5 is already included in the system- 
atic error covariance matrix, this has less than a 0.01 effect on 
w, about ten times smaller than it would have if we did not 
include this systematic. 

Next, we consider systematics caused by potentially differ- 
ent populations of supernovae. We p erform a cut on the best- 
fit true xi or c of each supernova (see lAmanullah et"an (120101) 
for details). The cosmology in each case is compatible with 
the cosmology derived from the whole sample. 

We now look at each of the four largest datasets for evi- 
dence of tensi on. Th e only tension found is in the first-year 
SNLS sample (lAstier et al. 20(51). Here, /3 and 5 are both at 
odds with the whole sample, but as noted above, we do not 
believe this is a cause for concern. 

The final two rows show th e hig h-redshift sample split by 
host type; this is discussed in 



5. CONSTRAINTS ON DARK ENERGY 

Following lAmanullah et al.l (120101) . we constrain the prop- 
erties of dark energy first using SNe la alone (with and with- 
out systematics), and then by combining the constraints de- 
rived from SNe la with those derived fro m the 7-year WMAP 
data of the CMB ( iKomatsu et alj|2011h . the position of the 
BAO peak fro m the combined anal ysis of the SDSS DR7 and 
2dFGRS data dPercival et al.ll20l"ol) . and t he measurement of 
the Hubble constant {Hq) from Cepheids (iRiess et al.ll201 lb . 

The rate of expansion at redshift z, H{z), is described by 



the Friedman equation: 



Hi 



(6) 



Hue cxp 



3(l + io(z))dln(l + , 



where Hq is the rate of expansion today, n„i and Hue are 
the matter and dark energy density with respect to the criti- 
cal density today, u'(z), is the dark energy equation-of-state 
parameter, and ilk = 1 — — ^ue is the spatial curva- 
ture density of the universe. Distances, such as the luminosity 
distance, depend on the integral of 1/H{z) over redshift. 

In this section, we consider the following models for dark 
energy; 

- ACDM: A cosmological constant in a flat universe. 

- wCDM: A constant equation-of-state parameter in a flat 
universe. 

- owCDM: A constant equation-of-state parameter in a 
curved universe. 

- w^CDM models: A time-varying equation-of-state pa- 
rameter in universes with and without curvature. 

The results for each of the models are listed in Table |7] 
and discussed in turn in the following sub-sections. Unless 
stated otherwise, the uncertainties represent the 68% confi- 
dence limits (Ax^ ~ 1) and include both statistical uncer- 
tainties and systematic errors. 

5.1. ACDM 

In the ACDM model, the equation-of-state parameter is ex- 
actly — 1 and does not vary with time. In a flat Universe, 
SNe la alone constrain the dark-energy density, Ha, to be 



0.705 



+0.040 
-0.043 



including systematics. Adding the con- 



straints from CMB, BAO and Hq reduces the uncertainty. Un- 
der the assumption of a flat Universe, the four probes yield 



Ha = 0.729] 



0.014 



(ACDM : SN + CMB + BAO + Hq) 



In this ACDM model, the expansion of the universe 
switched from deceleration to acceleration at z = 0.752 ± 
0.041, which corresponds to a look back time of 6.62 ± 
0.22 Gyr, about the half of the age of the universe. Equal- 
ity between the energy density of dark energy and matter oc- 
curred later, at z = 0.391 ± 0.033 or the look back time of 
4.21 ±0.27 Gyr 

In Figure |5] we show the confidence intervals on and 
Qa from SNe, CMB and BAO. Both the individual constraints 
and the combined constraint are shown (the BAO constraints 
are computed with an il,nh^ prior from the CMB). The SN 
constraint is almost orthogonal to that of the CMB. If we re- 
move the flatness prior, the best-fit ilm and Ha change by a 
fraction of their errors with Qi. ~ 0.002 



+0.005 
-0.005- 



5.2. wCDM : Constant Equation of State with to 7^ —1 

In wCDM models, w is constant but is allowed to be differ- 
ent from —1. W hile few dark energy theories give w ^ —1 
and yet constant ( ICopeland et al.ll2006l) . the constant w model 
is still useful to constrain as it contains fewer parameters than 
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Table 6 

Constraints on standardization and cosmological parameters for subsets. Mg is the _B-band corrected absolute magnitude; a, /3, and 5 are 
the lightcurve shape, color, and host mass coiTection coefficients, respectively. The outHer rejection is redone each time, so the totals may not 
add up to the whole sample. The constraints are computed including BAO, CMB, and Hq constraints and supernova systematic eiTors. 



Subset Number MbQi = 0.7) a fi 5 U.m w 



Whole Sample 



z > 0.015 


580 


-19 321+0 030 


121+0-007 
o. J-ZJ-_0.007 


2 47+0.06 


-0 032+0-031 

U.UOZ q 031 


r, 97-1+0.015 
'J-^' -1-0.014 


1 01-^+0-068 
l.uio_o 074 


Correction Coefficients, Split by Redshift 


0.015 < 2 < 0.10 
0.100 < z < 0.25 
0.250 < z < 0.50 
0.500 <z< 1.00 
z > 1.000 


175 
75 
152 
137 
25 


-19.328lO;°l 
-19.37linK 

-19 S17+0 046 
i9.JU(_Q 049 

-19.289l°;il 


n 11S+0-011 
U.il8_o.oii 

n 146+0-019 

U.i4D_o 019 

n 116+0014 

U.ilD_o 013 

n 124+0-019 

U.iZ4_o 019 
-0 019+0 072 
U-Uiy_o 076 


2 C.7+O.O8 
-0.08 

2.56l°;i? 
3-48_o;89 


-n 027+0-05-1 
'J-U^' -0.054 
n 087+0-060 

-U.0»r^_0,o60 

-n 042+0-066 
U.U4Z_0 066 
023+0-060 
U.UZJ_o 060 

-n 151+0-384 
U.101_o 446 


0.270 (fixed) 
0.270 (fixed) 
0.270 (fixed) 
0.270 (fixed) 
0.270 (fixed) 


-1.000 (fixed) 
-1.000 (fixed) 
-1.000 (fixed) 
-1.000 (fixed) 
-1.000 (fixed) 








Effect of 5 


on w 








z > 0.015 
z > 0.015 


580 
580 


-19 s4n+o o26 

iy.J'iU_Q 026 

-19 sn3+°-o^^ 


r. ,23+0-007 
U.1Z,3_0 007 
120+0-007 


2 47+0.06 
^-^'-0.06 
2 47+0.06 


-0.080 (fixed) 
0.000 (fixed) 


272+0.015 
U.Z(Z_o 014 
n 97-1+0.015 
'J-^'l-0.014 


-1 004+0-067 
i.UU4:_o 072 

, niQ+o-oeg 

l.U10_Q 075 



Cosmological Results, Split by Lightcurve Color and Shape 

'=>0-05 ii^ -I9.387tlill O.llSlO;"!! 2.77tro^9 -0-057lO:gg 0.2691^1! -1.028«:«^I 

c< 0.05 321 -19.323to«^3« 0.125to«lS 1.29l«:g -0.057«:°i 0.2751^1! -0.982to.069 

XI > -0.25 311 -19.366lO-«« OmO+^ilt 2.58+_'-]° -0.004^-047 o.269+_°-?}A -1.037t° ?J' 



-0.041 U.UZU_o.o25 ^-^^-0.10 -^•>J^'t-0.047 ^-^O»-0.014 

XI < -0.25 269 -19.386tO;0*^ 0.152l0-02i 2.43l°:°8 -0.087tO;«^o 0.267t^,f,l -1.045^ 



Con'ection Coefficients and Mg for the Large Datasets 



Hicken et al. (2009) 94 -19-314^0 055 0.115to oi5 2.74toH -0.053lo o9| 0.270 (fixed) -1.000 (fixed) 

Holtzman et al. (2009) 129 -19.3361 05I O.Ugio olg 2.40lo;i4 -0.06lto.oto 0.270 (fixed) -1.000 (fixed) 
Miknaitis et al. (2007) 74 -19.325innRO 0.113lno35 2.49ln!fi 0.000 (fixed) 0.270 (fixed) -1.000 (fixed) 



''-O.OSO '-'-ll^''-0.035 ^-^='-0.16 
Astier et al. (2006) 71 -19.292l^-!54g 0.145tool8 l-70loig -0.023l:oolo 0.270 (fixed) -l.OOO(fixed) 



z > 0.9, Split by Galaxy Host 

Early Type 2 > 0.9 13 -19.388^ oilse ^■'^^'^-O.ISI ^■'^^-l.lt 0.000 (fixed) 0.270 (fixed) -1.000 (fixed) 

Late Type 2 > 0.9 15 -Id.Ulto'^fr 0.094l[j;5340 0.49lo;g| 0.000 (fixed) 0.270 (fixed) -1.000 (fixed) 



the dynamical dark energy models considered in the next sec- 
tion, and a value different from — 1 would rule out the cosmo- 
logical constant. 

In a flat universe (ilfe= 0), SNe la alone give w = 
— l.OOl^Q ^gg (including systematics). Adding the constraints 
from the other three probes tightens the constraint on w con- 
siderably, as the constraints from SNe la in the ilm- w param- 
eter plane are almost orthogonal to those provided by BAO 
and the CMB (Figure |6]l. 

In principle, a constraint on Hq helps to break the degen- 
eracy between and h for CMB, which measures Vtmh? 
jSpergel et al.ll2003l) . However, in this case adding supernova 
data helps more, as narrowing the degeneracy between ilm 
and w allows the CMB itself to constrain Hq. By combin- 
ing all four probes, we find w ~ — 1.013j^Q Q73. As seen in 
Table |7] neither BAO nor i/o currently make much of a dif- 
ference in the error bars for this model. 

5.3. owCDM : Constant Equation of State in a Curved 
Universe 

Inflation models general ly predict that the curvature of th e 
Universe, flk, is - IQ-^ (iGuthll 19811: Liddle & LvthlllOOOh . 
In curved universes, SNe la play the critical role in con- 
straining w, while CMBh-BAO constrain ftk and flm- By 



combining all four probes, we find Qk = 0m2Za"taj and 
w = — l.OOSj^Q Qgg. Even with the additional freedom for 
non-zero curvature, a flat universe is supported from observa- 
tions. Among many cosmological parameters, the curvature 
of the universe is the most well-determined parameter 

We note CMB alone does not place a tight constraint on 
curvature', Vlk = -O.l02ta"olj (iKomatsu et alj|2011h . and 
it is the combination of SN la and BAO which improves the 
constraint by a factor of ten. Including Hq improves the cur- 
vature constraints but the current measurements from SNe la 
and BAO have more impact. 

For the equation-of-state parameter w, as seen in Table |7] 
BAO constraints are now needed to constrain 0„i, while Hq 
again has a small impact on the w measurement. 

5.4. Time Dependent Equation of State 

We next examine models in which dark energy changes 
with time. For a wide range of dark energy models, it can 
be shown ( lLinded2003h that, to good approximation, the dark 
energy equation-of-state can be parametrized by 

w(a) ^ Wo + Wail - a) (7) 

* |http : //lambda ■ gsf c ■ nasa ■ gov/product/map/dr4/parameters ■ cfiii| 
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Table 7 

Fit results on cosmological parameters Qm, wg, Wa and $1^. The parameter values are followed by their statistical (first column) and statistical and systematic (second 
column) 1(7 (Ax^ = 1) uncertainties. For the fits including curvature and time-varying w, the confidence intervals can be quite non-gaussian and we also show Ax^ = 4 

confidence intervals (with and without systematics) for comparison. 



Fit 




Qm w/ Sys 


fife 


Qfe w/ Sys 


wo 


•(OO w/ Sys 




Wa w/ Sys 




ACDM 


SNe 


r. 977+0.022 
'-0.021 


U.^»0_Q Q4Q 


(fixed) 


(fixed) 


-1 (fixed) 


-1 (fixed) 


(fixed) 


(fixed) 


SNe+BAO+CMB 


n 97S+0-014 
U-^'»_0.013 


n 9H9+0.017 

U.ZOZ_Q gj^g 


(fixed) 


(fixed) 


-1 (fixed) 


-1 (fixed) 


(fixed) 


(fixed) 


SNe+BAO+CMB+Ho 


^•'^'^-0.012 


n 271+°"" 
U-^'l_0.014 


(fixed) 


(fixed) 


-1 (fixed) 


-1 (fixed) 


(fixed) 


(fixed) 


oACDM 


SNe+BAO+CMB 
SNe+BAO+CMB+Ho 


r, r,on+0.015 

r, 971+0.013 
'^■^'^-0.012 


0.286lO:°l? 

r, 979+0.014 
"■^'^-0.014 


-0 004+°-'""^ 




-1 (fixed) 
-1 (fixed) 


-1 (fixed) 
-1 (fixed) 


(fixed) 
(fixed) 


(fixed) 
(fixed) 


wCDM 


SNe 
SNe+BAO+Ho 
SNe+CMB 


n 901 +0.067 
U.JUy_Q 

r, 971 +0.018 


n •?2n+°-'^35 

U.OZU_Q 

0.279lO;°r3 


(fixed) 
(fixed) 
(fixed) 


(fixed) 
(fixed) 
(fixed) 


i.Uli_Q 231 


i.UUl_Q 3gg 

i.U(D_Q 433 
-n 9';^+0 075 
U.»00_o.o79 


(fixed) 
(fixed) 
(fixed) 


(fixed) 
(fixed) 
(fixed) 


SNe+BAO+CMB 
SNe+BAO+CMB+i?o 


n 97S+0-014 

r, 979+0.013 
"■^'^-0.013 


0.285ini? 
0.27linil 


(fixed) 
(fixed) 


(fixed) 
(fixed) 


U.M»J_o.055 
-i.UU8_o 054 


-n 9';i+0-0''5 
u.yoi_o.o8i 
1 niq+0.068 

±.U10_Q Q73 


(fixed) 
(fixed) 


(fixed) 
(fixed) 


owCDM 


SNe+CMB 
SNe+BAO+CMB 
SNe+BAO+CMB+i?o 


r. 901 +0.069 
U-^oJ^-0.087 
n 9sq+0-016 

r. 979+0.013 
"■^'^-0.013 


0.295lO;l°? 
0.287in}? 

r, 979+0.015 

U.Z(Z_Q 


U.UUO_Q Q27 


-n 002+" °°* 
n nn2+°°'^'' 

U.UUZ_g QQ7 


l.UU(_Q ^94 
1 ni9+0'058 

1.U1Z_Q Qg2 
i.UUD_Q QgQ 


-0.993t«jl? 
-0.975t-.E!r8 


(fixed) 
(fixed) 
(fixed) 


(fixed) 
(fixed) 
(fixed) 


WzCDM 


SNe+CMB 
SNe+BAO+CMB 
SNe+BAO+CMB+J?o 


r, 97,,+0.022 

U.Z IO_Q Q20 

n 97S+0.014 

U-^'°-a.oi4 
n 97-1+0.013 
'J-^'^-o.ois 


n 901+0.043 
u.zoi_Q 028 

0.284t0;«l? 
0.270l0;«il 


(fixed) 
(fixed) 
(fixed) 


(fixed) 
(fixed) 
(fixed) 


1 nnc+0.165 

-l-006_o.i82 
-1 n52+"-i26 

-1 021+0-123 
l.UZi_(j 


-n 9q'^+"-263 
U.jyj_Q 307 

1 rti Q+0. 183 
l.UiO_Q 473 

i.U4D_Q 470 


„ -,1+0.75 
'J- -0.77 

'^•'^^-0.62 

n 07+0-49 

-0.60 


n 17+1-08 

'J-i'-i.ig 
n 26+0-5'^ 

U-ZD_o 74 

n 1 4+0-60 

U-14_o 76 



ow-CDM 



n 07^+0.065 
U.U/0_Q 42g 

-n 004+0-017 
U.UU4_Q 010 

02^+0-008 
U.UZO_Q 008 

02'5+0-016 
U.UZO_Q^035 



0.073t°-i4if 
-0 001+0-037 

U.UUi_o 013 

027+0-012 
U-UZ/_Q oil 

027+0-026 
U.UZ/_Q^036 



SNe+CMB 
SNe+BAO+CMB 
SNe+BAO+CMB+ffo 
SNe+BAO+CMB+,f/o (Ax^ = 



4.0) 



n 1 77+0.086 

"■1 ' ' -0.093 

n 9Q-3+0.019 

U.^83_0.017 

270+0-014 
U-Z(U_o 013 

270+0-029 
^•^'*J-0.026 



0.190iO-.?»« 



0.286lO-«f3 

^+0.016 
-0.015 



0.274"' 



rs 971+0.032 

^•^'*-0.029 



0.17C 
-0.202 
0.169 
0.178 
0.069 
0.072 
-, 9-10+0.425 
-l.^l»_0.147 



-0.988 
-1.010 
-1.218 



-0.969"' 



■1+0.284 
-0.345 
-0 997+0-266 

u.yy/_o 293 

1 -IQO+O.IOO 

-i.iys_o 442 

1 -igo+0.293 
-i.iy«_0.227 



0.90 
-0.01 
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where a = l/(l-|-^;)isa scale factor. The ACDM model is 
recovered when wq = —1 and Wa = 0. The constraints on wq 
and Wa are shown in Figure |7] and Table [7] 

The Dark Energy Task Force (lAlbrecht et al.l |2006|) pro- 
posed a figure of merit (FoM) for cosmological measurements 
equal to the inverse of the area of the 95% confidence contour 
in the wq — Wa plane. When we make this measurement, using 
the Ax^ — 5.99 contour, we find a FoM of 1.84 (statistical- 
only) and 1.04 (including systematics). Frequently, the FoM 
is also defined in terms of the la errors (Ax^ = 1); this 
FoM is 39.3 (statistical-only) and 22.6 (including systemat- 
ics). Surprisingly, even with Wa floating, we still find an in- 
teresting constraint on of ~ 0.02. 

We next consider a model in which the dark-energy 
equation-of-state parameter is constant inside fixed redshift 
bins. This model has more parameters (and thus more free- 
dom) than Wq — Wa- The results are shown in Figured and 
Table [7] We adopt the redshift bins used in Amanullah et al. 
(1201 Oh . so that a direct comparison can be made. 

In the left panel with broad bins, we show a reasonably 
good measurement of the equation-of-state parameter from 
redshift to 0.5. From redshift 0.5 to 1, there is no real 
constraint. For example, any scalar field model {\w\ < 1) 



is reasonably compatible with the data. Above redshift 1, the 
constraints are weaker w > is ruled out, as this violates 
early matter domination. 

We separate the supernova and early universe constraints by 
defining a bin at redshift 1.6, as shown in the middle panel. 
This shifts the confidence interval for w(l.G < z < 1.6) to- 
wards higher w. Eliminating this division, and instead adding 
more bins up to redshift 0.5 (right panel), gives three con- 
straints of moderate quality with a possible crossing of w = 
— 1. No matter the binning, we will need more data extending 
above redshift 1 to investigate the dark energy equation-of- 
state parameter where the uncertainty is still very large. 

To examine constraints on the existence of dark energy at 
different epochs, we study p{z), which is the density of the 
dark energy and allowed to have different values in fixed red- 
shift bins. Within each bin, p is constant. (Note that the dis- 
continuities in p{z) at the bin boundaries introduce disconti- 
nuities in H{z).) We choose the same binning as above, but 
note that binned p and binned w models give different expan- 
sion histories. Our results are shown in Figure 0and Table 

m 

Although there is no real constraint on the equation-of-state 
parameter at redshift 0.5 to 1, dark energy is seen at high sig- 
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Figure 5. ACDM model: 68.3%, 95.4%, and 99.7% confidence regions of the {Qm, SIa) plane from SNe la combined witli the constraints from BAO and 
C]\4B. The left panel shows the SN la confidence region only including statistical errors while the right panel shows the SN la confidence region with both 
statistical and systematic errors. 




Figure 6. luCDM model: 68.3%, 95.4%, and 99.7% confidence regions in the (Q™, «)) plane from SNe la BAO and CMB. The left panel shows the SN la 
confidence region for statistical uncertainties only, while the right panel shows the confidence region including both statistical and systematic uncertainties. We 
note that CMB and SN la constraints are orthogonal, making this combination of cosmological probes very powerful for investigating the nature of dark energy. 
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nificance in both panels. There is weak evidence for the exis- 
tence of dark energy above redshift 1, as can be seen in the left 
panel. However, if we again separate the supernova data and 
early universe constraints (right panel) we see neither probe 
has any constraint on the existence of dark energy above red- 
shift 1. 

6. DISCUSSION 

6.1. Improving the Constraints on Time-Varying w by 

Efficiently Adding z > 1 Supernovae 

Beyond z = 1, we add 10 new well-measured SNe la to 
the Hubble diagram. The variance-weighted RMS scatter of 
the new sample is 0.20 ± 0.05 mag. As a comparison, the 
15 z > 1 SNe la from the GOODS survey that pass our 
Union2 selection cuts have a variance-weighted RMS scat- 
ter of 0.25 ± 0.05 mag. The new sample almost doubles the 
weight of HST-discovered SNe la beyond z = 1. The increase 
provides improvements on the most difficult-to-measure pa- 
rameters, those that describe the time-varying properties of 
dark energy: p{z) and 'w{z) at the higher redshifts. In particu- 
lar, the supernovae from this search improve the constraint on 
p{z) at redshifts l.Qz < 1.6 by 28% (statistical errors only) 
and 18% (including supernova systematics) after adding the 
constraints from the CMB, BAO and Hq (using the binning 
illustrated in the right panel of Figure|9}. (It is more difficult 
to compare binned w results, as the constraints are much less 
gaussian and more sensitive to the location of the best fit.) 

The new sample is also obtained with greater observing ef- 
ficiency with HST. Considering the number of z > 1 SNe la 
that make the Union2 selection cuts, the yield of SNe la in- 
creases from a rate of one SN la per 43 HST orbits in the 
GOODS survey to one SN la per 22 HST orbits in this survey. 

6.2. Splitting the sample according to host galaxy type 

SNe la are well-standardized with a small dispersion in 
magnitudes across the whole class. Any clues to heteroge- 
neous characteristics therefore offer exciting possibilities to 
further improve standardization, enhancing the use of SNe as 
a cosmological probe. There is now evidence from studies 
of large samples of SNe la at both low and intermediate red- 
shifts (0 < z < 0.8) that SN la properties are related to the 
properties of the host. The clearest of these is the relation be- 
tween light curve width and the specific star formation rate. 
SNe la in passive galaxies tend to have narrower light curves 
than SNe la that are in galaxies that are actively forming stars. 

More than two-thirds of our new SNe la beyond z = 0.9 
are hosted by early-type galaxies (Meyers et al. (2011)). In 
field surveys, such as the GOODS survey, this ratio is in- 
verted. By combining SNe la from our HST Cluster SN 
Survey and GOODS, t ogether with our z > 0.9 SNe la in 
lAmanuUah et al.l ( 120101) . which have HST images of the host, 
we can create a sample of SNe la that has roughly equal num- 
bers when split according to host type. When split this way, 
we find that z > 1 SNe la in early galaxies rise and fall more 
quickly than SNe la in later host types, thus extending the red- 
shift interval over which the effect is now detected. Finding 
that low and high redshift SNe la follow similar trends gives 
us confidence that we can use very distant events to constrain 
cosm ological pararneters. This finding is reported in more de- 
tail in lMevers et all (1201 lb . 

There is also evidence from SNe la at low an d intermediate 
redshi fts for other correlations with host type. ISullivan et al.l 
(120 lOi) find that both /3 and the RMS scatter about the Hubble 



diagram are smallest for SNe la in passive galaxies. These 
trends suggest that dust plays a greater role in reddening and 
dimming SNe la in late-type galaxies. We examined our z > 
0.9 sample for evidence of similar correlations using our host 
classification from Tables 3 and 4 of Meyers et al (201 1). 

After correcting SN la luminosities for lightcurve shape, 
SN la color and host galaxy mass (with the global values of 
these correction coefficients), we measure a sample disper- 
sion of 0.14j;Q |!)g mag for SNe la in early-type galaxies and 
0.14tQ Qg mag for SNe la in late-type galaxies. In terms of 
the RMS, we find 0.23 ± 0.05 mag and 0.26 ± 0.05 mag for 
early and late-type samples, respectively. The uncertainties 
are currently too large to distinguish between the two sam- 
ples. Similarly for (3, the errors are larger than the difference 
between the two samples, as seen in Table |6] Clearly, higher 
quality data of a larger number of z > 0.9 SNe la in both early 
and late-type galaxies are required before the trends that are 
seen at low redshift can be detected in high redshift samples. 

We also examined the error-weighted difference in the 
brightness of SNe la in the two samples after correcting for 
lightcurve shape and color, but without correcting for the host- 
mass luminosity relation (setting 5 = 0) and find that SNe la 
in early-type galaxies are 0.18 ± 0.09 mag brighter Since 
early-type galaxies are typically more massive than late-type 
galaxies, this 2a difference, if confirmed with larger statistics, 
may be related to host galaxy mass. 

6.3. Future directions with current instrumentation. 

Due to the much improved sensitivity of the WFC3 IR de- 
tector, it will be feasible to measure z > 1 SNe with much 
better precision. The color measurement errors 0.03 
in B — V") can be made comparable to t he color measure- 
ment errors in the SDSS supernova survey (iSmith et al.1120021 
iHoltzman et al. 2008). Assuming that the intrinsic dispersion 
of SN la luminosities does not change with redshift, the vari- 
ance weighted RMS of the WFC3 sample should be similar 
to that measured for the SDSS, i.e. ^ 0.14 mag. A well- 
observed SN la with WFC3 should have a statistical weight 
of two to three SNe la from the Cluster and GOODS surveys. 

With a sufficient number of well-measured z > 1 SNe la 
with WFC3, it should be possible to search for the correlations 
between the properties of SNe la and their hosts that are seen 
at lower redshifts. As discussed above, current samples at z > 
1 are too small to detect these differences. With the improved 
WFC3 photometry, only 40 SNe la, split evenly between early 
and late-type hosts, would be needed to constrain a difference 
in /3 to an uncertainty of 0.4, whi ch is about half the di fference 
found for lower redshift SNe la (ISullivan et al.ll2010l) . These 
samples would be just enoug h to see evidence of the lower 
RMS for passive hosts seen bv lSumvanetaLlJIoTol) . 

Current WFC3 SN la surveys target empty fields, which 
means that there will be few SNe la in passive host galaxies. 
A WFC3 SN la survey that spends part of its time targeting 
z > 1 clusters would ensure a better balance between host 
types while increasing the overall yield. 

In order to investigate the figure of merit constraints pos- 
sible with WFC3, we simulate a sample of 40 supernovae at 
redshift 1 .2 and add this sample into the current compilation. 
As there is a hard wall at + Wa = when including BAO 
and CMB data, we simply fix fi™, rather than including BAO 
and CMB data (the alternative would be to adjust the super- 
nova magnitudes to a cosmology model far away from the 
wall). When adding these supernovae, the statistical figure 
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Table 8 

Constraints on redshift binned equation of state w and density p (normalized by the current 
critical density). The consfraints are computed including SNe, B AO, Ho , and CMB data. 







z < 0.5 


0.5 < z < 1.0 


1.0 < z < 1.6 


z > 1.6" 


w(z) 


Stat Only: 
w/ Sys: 


l.U10_,j Qgg 

-1 006+""" 
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Note. — Constraints on binned pdeC^) / PcO ™d ui(2). This redshift binning con'esponds 
to the middle panel of Figure[8]and the right panel of Figure|9] 

" We note that the weak constraints in these bins come mostly from the CMB (which tells us 
that the early universe was matter-dominated) and are only indirectly constrained by super- 



of merit improves by 39%. By the same metric, the current 
cluster sample improves the figure of merit by 10%. 

6.4. Reducing the Systematic Errors for Future Surveys 

As has been stressed by several authors, systematic errors 
are now larger than statistical errors. To fully utilize the poten- 
tial of current and future SN la surveys to constrain cosmol- 
ogy, it will be necessary to reduce these errors significantly. 

The largest current source of systematic uncertainty is cal- 
ibration. Calibration uncertainties can be split into uncer- 
tainties related to the primary standard, and uncertainties in 
instrumental zeropoints and band passes. In principle, all 
of these uncertainties can be reduced by establishing a net- 
work of well-calibrated st andard stars and mon itoring tele- 
scope system throughputs (iRegnault et al.ll2009l) . The Sloan 
Digital Sky Survey demonstrated that a 1% relative photo- 
metric calibration is possible with the curr ent standard sta r 
network and system throughput monitoring (iDoi et al.ll2010h . 

The ongoing Near by Supernova Factory (SNf) project 
jAldering et al 1 120021) is aiming to provide the network of 
standard stars. SNf will also address the systematic uncer- 
tainty due to host-mass correction since the range of host 
masses would become comparable to that of high redshift for 
the first time. Additionally, the comprehensive SN la spectral 
time series from the SNf will allow one to tackle systematic 
uncertainties related to modeling of the lightcurves. 

In the future, recently approved experiments such as AC- 
CESS (Absolut^ColorCalibration Experiment for Standard 
Stars iKaiser et all 120 1 Ol) and the proposed NIST STARS 
project (National In st itute for Standards and Technology 
livicGraw et ail I20T0I: IZimmer et all I20T0I) are aiming to 
achieve sub-percent absolute flux calibration for the network 
of stars in the wavelength range of visible to NIR. With this 
network of stars and with new techniques for monitoring 
throughput of the telescopes (IStubbs et al.ll2007h . we will be 
able to cross-calibrate systems and reduce the systematic er- 
rors below the statistical errors. 

7. SUMMARY AND CONCLUSIONS 

We present HST ACS, HST NICMOS and Keck AO- 
assisted photometry of 20 SNe la in the redshift range 0.63 < 
z < 1.42. The SNe la were discovered in the HST Cluster 
Supernova Survey, a survey run by the Supernova Cosmology 
Project to searc h for SNe la in field s centered on 25 distant 
galaxy clusters ( iDawson et al.ll200"9l) . 

We implement new techniques to improve the accuracy of 
HST photometry. In particular, for data taken with NIC- 
MOS, which samples the rest-frame B and y-bands of z > 1 



SNe la, we use a more direct, more accurate measure of the 
NICMOS zeropoint (Ripoche et al. 2011), and we remove 
the residual background that persists after standard process- 
ing of NICMOS data with the CALNICA pipeHne (Hsiao et 
al. 2010). For data taken with ACS WFC in the z^so fil- 
ter, we incorporate a SED-dependent aperture correction (see 
Appendix A). 

F ollowing the p rocedures outlined in lKowalski et alJ (120081) 
an d iAmamillahet al. (2010), we add our SNe la to the Union2 
compilation. Fourteen of the 20 SNe la of our supernovae 
pass the Union2 selection cuts. Ten of them are at z > 1. 
The strategy of targeting high-redshift galaxy clusters results 
in factor of two improvement in the yield per HST orbit of 
well-measured SNe la beyond z ~ 1 and a factor of three 
to five improvement for SNe hosted by early-type galaxies. 
For WFC3, with its smaller field of view, the advantage of a 
cluster search is even greater 

We use the new Union2.1 sample to constrain the proper- 
ties of dark energy. SNe la alone constrains the existence 
of dark energy to very high significance. After adding con- 
straints from the CMB, BAO, and Hq measurements, we pro- 
vide the tightest limits yet on the evolution of dark energy 
with time: Wa = 0.141q ^g. Our sample improves the con- 
straints on binned p by 18% in the difficult-to-measure high 
redshift bin, 1.0 < z < 1.6. Even with a time- varying WQ-Wa 
model, the universe is constrained to be flat with an accuracy 
of2%inl7fc. 

The results from this new cluster-hosted supernova sample 
point the way to the next steps that are now possible with the 
WFC3 on HST, an instrument that can obtain high signal-to- 
noise, multifilter SN la lightcurves at z > 1. The cluster 
approach, used in this paper, would make it feasible to build 
a significantly larger sample at these highest redshifts, evenly 
balanced between early and late-type hosts. With such a sam- 
ple, we can mitigate the effects of dust and evolution that may 
ultimately limit constraints on time-varying w. 
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Figure 7. 68.3%, 95.4%, and 99.7% confidence regions of the {wg, Wa) plane from SNe combined witli the constraints from BAO, CMB, and Hq, both with 
(solid contours) and without (shaded contours) systematic errors. Zero curvature has been assumed. Points above the dotted line (wo + Wa > 0) violate early 
matter domination and are disfavored by the data. 
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Figure 8. Constraints on w(z), where w(z) is assumed to be constant in each redshift bin, are plotted at the 68% probability level (Ax^ = 1). Each panel 
shows different redshift binning. The results were obtained assuming a flat universe for the joint data set of SNe, BAO, CMB, and Hq, with (dark/orange) and 
without (light/yellow) SN systematics. The middle panel takes a closer look at the 2 > 1 constraints, while the right panel shows the effects of w binning at low 
redshift. In this panel the best fit values of w cross w = —1 twice at low redshift, an unusual feature in dark energy models. We note that the ACDM model is 
consistent with our w{z) constraints for each of these binnings. 
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Figure 9. Redshift evolution of dark energy density: Constraints on p(z) are shown as a function of redshift, where p(z) is the density of the dark energy at a 
given redshift bin and assumed to be constant within the redshift bin. p(z) is normahzed by the critical density today (pco) and is plotted at the 68% probability 
level (Ax^ = 1). The results were obtained assuming a fiat Universe for the joint data set of SNe la, BAO, CMB, and Hq, with (dark/orange) and without 
(light/yellow) SN systematics. The two panels demonstrate different redshifts binning and have different scales. 
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APPENDIX 

APPENDIX : ACS COLOR DEPENDENT APERTURE CORRECTION 

The scattering of long wavelength photons (> 8000 A) within the ACS CCDs causes the point sp read function (PSF) of images 
taken in the ACS zgso-band to depend on the spectral energy distribution (SED) of the source (ISirianni et al.l 120051 iJee et al.l 
l2007h . The SED-dependent PSF means that aperture corrections also depend on the SED. In this section, we describe how we 
derive aperture corrections for our observations taken with the ACS Z8,5o-band. 

We use two stars (GRW70 and GD72 from ACS programs PID9020 and PID 10720) to measure the aperture correction (AC) 
as a function of wavelength. The stars were observed with 15 narrow-band filters between 7660 A and I0360A. We processed the 
data in the same way as the SN la images. To derive the aperture correction, we compare the flux in two apertures, one with a 
3-pixel radius and the other with a 1 10-pixel radius (the radius used for defining the zeropoint). Errors from removing the sky are 
the dominant source of the uncertainty of this measurement. 

The measured encircled energy and its best-fit curve for a 3-pixel radius are shown as a function of wavelength in Figure [Tol 
We apply an aperture correction to the observed flux with this best-fit curve. Figure [TT] demonstrates that the ^775- zsso color 
is not accurate enough for SNe la to infer the accurate aperture correction. Therefore, we need to introduce an SED-weighted 
aperture correction. 
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Figure 10. Wavelength vs. Encircled Energy (EE) Fraction for a 3-pixel aperture. The solid line is a fit (Eg. lAlt to the measured data points. The dotted line is 
the relation from Sirianni et al. (2005) 



We find that the best-fit encircled energy curve can be described with 



EE(A) 



1 



L28 + 1.45 X 10-8 X exp(17.99 x A) 



(Al) 



where A is wavelength in /im. The formula is valid in wavelength for virtually all the F850LP filter 

Since the SED of a SN la is a function of redshift and phase, so is the aperture correction. We demonstrate two methods to 
account for this dependence. The first method, "Method I," calculates the aperture correction iteratively using the SED produced 
by SALT2. The advantage of this method is that the final aperture-corrected fluxes can be used with the standard throughputs 
and zeropoints. The second method, "Method II," computes a new zeropoint and F850LP throughput for direct use with the flux 
measurements that have not been aperture-corrected. We note here that an apparent color difference between SNe with only ACS 
data and those that also had NICMOS observations was seen unblinded before the development of these methods, so the relative 
colors of the ACS-only and NICMOS -included subsamples should not be considered bhnded. 
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Figure 11. Wavelength vs. 4775 —2850 color: solid asteiisks are stars from Gunn-Strvke r Catalog ^Gunn & Strvkeil 19831) : open asterisks are stars from Pickles 
catalog iPickles 1998). The dotted hne is the best-fit line to the stellar locus. SN la loci jHsiao et al.ll2007l) are plot for three redshifts with color-coded phases 
and arrows which show the direction of the phase evolution. The SN la loci deviate from the stellar locus as the SED is different. Therefore, 4775 —^850 colors 
are not accurate enough to perform apeiture connections, necessitating the introduction of an SED- weighted aperture coirection. 

Both methods give the same answer, ahhough the error in the color from "Method 11" is sHghtly larger than that of "Method I", 
because the effective wavelength of the 2350 shifts towards the 1775 filter, thus shortening the wavelength separation of the two 
filters. For this paper, we adopt "Method 11" and report all results using this method. 

Method I: Iterative approach 
The magnitude of an object in the Vega magnitude system is defined (iFukugita et al.|[T995t ISirianni et al.ll2005h as 

where A is wavelength, R is the system response, fx is the SED of the object (e.g. SN la), and /AVcga is the SED of Vega. We 
define the Vega magnitude zeropoint (Zpt) as 

Zptvoga = 2.5I0510 (^J XRfxvcgiidX^ . (A3) 

Eq. lA2l can then be rewritten as: 



m = -2.51ogio XRfxdXj + Zptvega- (A4) 

The observed magnitude, niobs, within a given aperture is 

mobs = -2.51o.gio (^J X REE{X) fxdX^ + Zptycga, (A5) 

where EE(A) is the wavelength dependent encircled energy (EE), which can be derived from Eq. lAll 
We define the SED dependent aperture correction (Atocoit) as 

f f XREE(X) fxdX\ 
= -2.51o,,o jjf^^i ) (A6) 

In practice, we do not know fx in advance, so we cannot compute the aperture correction directly. Instead, we derive it 
iteratively using the SED derived from fitting the SN la light curve with SALT2 as input to the next iteration. With this method, 
we use the STScI filter response function and zeropoints. 

Method II: Modified Filter with Zeropoint 

We can rewrite Eq. lA5l as: 



-2.51ogio (^jl^^^^^lM^^ + Zptvega - 2.5I0510 (^J XREE{X) fxve^.dX^ . (A7) 
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Effectively, the last term serves as a zeropoint offset for a given aperture radius. 
We rewrite Eq. lA7l as: 

/ [XREEiX)fxdX \ 
mobs = -2.51ogio r > t,t:t:/\\ c 7T + Zptvcga - AZptvoga, (A8) 



/Ai?EE(A)/AVcgarfA 

AZptvcga = 0.438 for a 3-pixel radius aperture and the zgsg filter We interpret Eq. |A8| as a magnitude measurement that uses 
a modified filter response, i? EE(A), and an adjusted zeropoint, Zptvcga — AZptvoga- We then run SALT2 using the counts in a 
fixed aperture, the modified filter response and the adjusted zeropoint. 
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